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ABSTRACT 


This thesis investigates the feasibility of adapting Naval 
Architecture and Marine Engineering design processes to digital computer 
solution. 

Two problems, the determination of critical whirling frequencies 
of propeller shafts oA the determination »f stresses in a transverse web 
frame are programmed for a digital computer, A detailed explanation of 
the programs and a discussion of the methods of approach and difficulties 
encountered in adapting ah solutions of these problems to a digital 


computer is included. The most important results obtained from programming 





these problems are experience with and knowledge of digital comm vers. 
Applying this knowledge and experience to the general field of naval 


architecture the following conclusions and recommendations are presented: 


Conclusions: 
a. Digital computers are applicable to naval architecture problems. 


b. The average naval architect has sufficient background to pro- 
gram a digital computer. 


ec. A programmer should have a comclete engineering understanding 
of the problem and some familiarity with numerical methods of 
mathematical processes. 


d. in general the limitations imposed by the computer on a 
problem are no more restrictive than those imposed by manual 
methods. 

e. The digital computer used should be able to handle numbers 
usually encountered in engineering problems, 

Recommendations: 

a. Digital computer should be used in naval architecture, 

b. Digital computers are suited for theses and research problems. 

c. Any future thesis should only consider the computer as an aid 
in the solution of a specific problem and not general appli- 
cability of the computer. 

ad. For engineering applications the programmer should have; (a) a 


good engineering understanding of his field; (b) knowledge of 


numerical methods of mathematical processes; (c) organizational 
abilities, 


Thesis Supervisor: Robert J, Hansen 


Title: Associate Professor of Structural Engineering 
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1. 


Introduction 





The field of naval architecture contains some of the most 
complex problems facing the engineering profession in this modern day 
and age, A ship is a structure containing a network of mechanical 
equipage that is unequaled by any other movable structure in the world. 
The complexity of such a structure almost defies analytic approach and 
as a result much of the design work is based on previously built vessels. 
A parent design is altered by repeating various numerical procedures 
until the desired result is obtained. The practice of naval architecture 
can best be described as laborious and repetitive with a necessity for 
a background of experience. 

In an effort to eliminate the necessity of utilizing personnel 
to perform laborious calculations, following World War II a committee was 
appointed by the Bureau of Ships to investigate and report on the appli- 
cability of utilizing machine aided calculations for problems arising in 
the preliminary design stage. Due to the stage of development of elec- 
tronic computers at that time an analog computer was developed for 
Stability calculations in the preliminary design. Although this investi- 
gation did consider the application of digital computers to problems in 
the preliminary design stage, the application to the general field of 
naval architecture was not fully considered, eho 

The purpose of this thesis is to investigate the applicability 
of digital computers to the field of naval architecture. with primary 


emphasis on structural design problems, to acquaint the authors and other 


1. Numbers in parentheses refer to the Bibliography at the end of the 
thesis. 
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naval architects with this new method of approach, to determine the type 
problems best suited for computer solution and publish some of the diffi- 
culties encountered, to determine the type of person required to utilize 
most effectively the machine aids and to incite interest in applying 
digital computers in the investigation of problems heretofore considered 
too laborious and time consuming to analyze by hand calculations, To 

do this the authors had to first acquaint themselves with an existing 
computer and secondly, to pick sample engineering problems of generic 
nature for specific programming. 

The authors are graduate students in naval architecture whose 
studies have been primarily channeled in the direction of hull structures. 
The background in mathematics included calculus through differential 
equations and a basic course in vector and complex mathematics, Their 
background in mathematics is no better than can be expected from the 
practicing naval architect, 

Two problems were picked as samples, one, calculation of the 
whirling frequencies of a propeller shafting arrangement and two, an 
analysis of an indeterminate structure of five degrees of indeterminancy. 
These problems were picked because they are examples of types that require 
reiterative procedures in their solutions. Further, both problems had 
solutions that had been previously determined by hand. 

The digital computer investigated and utilized was the Whirlwind I 
computer at MIT, This computer started operations in November 1950 and 
has been used as a guide for the development of many commercially available 
machines today. For this reason the general conclusions drawn from an 
analysis of the applicability of Whirlwind I should apply to any machine 


considered for particular problems, although other machines may have different 


capabilities. 
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Basically, all digital computers utilize the same theory of 
operation. They all contain the following elements: storage element, 
control element, arithmetic element and input and output devices, as 
depicted in Figure lI. 

Calculations are performed in the arithmetic element utilizing 
the numbers and instructions contained in the storage element. Input 
and output devices are used to get information in and out of the computer. 


The control element makes the components function. 


ARITHMETIC CONTROL SRORAGE 


ELEMENT x ELEMENT ELEMENT 


DIGIT TRANSMISSION BUS 





Schematic Diagram of Digital Computer 


Figure I 
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The discussion of the individual elements that follows pertains 
to Whirlwind I only, although other computers have similar equipment that 
performs the same functions, 

The storage element consists of 2048 registers or "cells." 
Instructions or numbers are stored in these registers by impressing 
electrical pulses on magnetic cores. Each register contains 16 such core 
and, therefore, the basic word-length is 16 digits. The pulses may be 
either positive or negative and the resulting direction of magnetization 
determines whether the digit represented is 0 or 1. This means that 
all the machine claculations are carried out in the binary number system. 
In every register the first binary digit represents the sign of the number 
and the remaining 15 binary digits can represent approximately 5 decimal 
lakes mics 

instructions are stored in the registers in a slightly different 
manner, Here the first five binary digits represent the mathematical 
operation that is to be performed when the instruction is obeyed. The 
remaining eleven binary digits represent the location of the register 
containing the number on which the operation is to be performed. 

The permanent storave of Whirlwind I may be increased to 
approximately 22,000 registers by the use of a magnetic drum, This device 
is somewhat slower in the transfer of information than is the permanent 
storage, but it is very convenient for the solution of long problems. 

The arithmetic element consists of three registers, In one of 
these, called the accumilator, or "AC" repister, the basic arithmetic 
operations of addition and subtraction are performed, Multiplication and 


division are accomplished by a series of these two basic operations. 
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The input devices are a photoelectric tape reader and a mechanical 
tape reader. These devices read into the storage element coded information 
from a punched paper tape. They convert holes in the tape into electrical 
signals to properly magnetize the registers of the storage element. A 
conversion routine is available to transform decimal information on the 
tape into binary information for storage in the computer. This feature 
is a great aid in programming and enhances the flexibility of the computer. 

The output devices available are flexowriters, five magnetic 
tape units and an oscilloscope and camera. The desired output information 
can be typed out directly by flexowriters, or it can be stored on magnetic 
tape and typed out by the flexowriter at a later time, The latter procedure 
Saves computer time since it is much faster to read information onto 
magnetic tape than to type it. Another method of retting the information 
out is to display it on the oscilloscope and photograph it. This saves 
typing time, but the film must be processed. However, it is convenient 
for reproducing the results. Standard output programs are available or 
special ones may be written to get the information out in the form desired. 
For instance, numbers may be in octal or decimal system; written in columns, 
lines, blocks, or tabulated; formed raised to a fixed power of ten or with 
a fixed decimal point, This flexibility of output is an important feature 


of digital computers. 


Control Element 


Rach of the various components of the computer consists of 
combinations of electrical circuitry and transfer buses. Only through 


intelligent control of these devices can electrical signals be made to 
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perform arithmetic operations. The control element performs this function 
by causing the flow of magnetic signals to reach the proper element at the 
appropriate time. Basically then, the control element is a switching 
device that channels the information through transfer buses to itself 

and the other elements. To do this, the control element must interpret 
the program instructions, which are themselves magnetic signals, into 
other magnetic signals which control the operation of the computer, 

The speed of operation of the computer depends upon the memory 
access time and the time required to integrate a number of switching com- 
binations into a useful operation, Switching is performed by utilizing 
electrical impulses from a high frequency oscillator. A new set of 
switching combinations can be established for each time pulse received 
by the control element. The time required to perform an operation depends 
on the number of impulses necessary to perform the switching and the 
repetition rate of the impulses, In Whirlwind I the pulse repetition rate 
is 1 micro=second. At this rate the addition of two numbers requires })8 
micro-seconds and multiplication requires aprroximately 65 micro~seconds. 
Using the conventional rating system which is the number of additions a 
computer can perform per second, the Whirlwind I is rated at 20,000 opera- 
bLONS per second, 

The instructions necessary to perform the desired mathematical 
calculations are called a program and are fully explained in Appendix A, 
For each computer there exists a standard instrvction code. By ingenious 
combinations of these standard instructions tre programmer can cause a 


comouter to perform the requisite numerical calculations, 
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One of the more important instructions is the "conditional 
response." This instruction allows the computer to make an either-or 
decision based on the sien of the number in the accumulator. Asa 
result of this capability the computer can be made to cycle back through 
a series of instructions any desired number of times, The cycling ability 
is ideally suited to repetitive processes which are useful in design calcv- 
lations, 

Digital computers have been widely used in other fields of 
engineering and new adaptations are being developed every day. Some of 


the more important uses developed are: 


(1) Aircraft weight, balance and thrust calculations, 
(2) Steam turbine design 
(3) Catapault launching analysis 


(4) General structural application in the aircraft 
industry. 





PROCEDURE 
General 


In order to accomplish the objectives discussed in the 
introduction it was necessary that some knowledge of computers be obtained. 
A twenty-five hour course in programming for Whirlwind I was completed 
by the authors to partially fulfill this requirement. 

Once the capabilities of the computer were realized, it was 
necessary to decide what problem would best prove the basic contentions, 
Probably the biggest single factor in this decision was time, A very 
complicated problem wovld require a great deal of time to program while 
a simpler problem would not test the authors' ability. 

A compromise was reached by deciding to program two problems. 
The first should be a relatively short and uncomplicated problem aimed at 
giving the authors experience with the computer and also making some 
contribution to naval architecture, The second problem should be of such 
Size and complexity that its solution would more fully describe the 
abilities of the computer and the benefits of computer applications to 
naval architecture. 

The first problem programmed was the determination of the 
whirling frequencies of propeller shafts. The method of solution was 
developed by Professor F, M. Lewis. This problem was selected because 
the method of solution is an iteration process ideally suited to computer 
application. Furthermore, the problem was short and each step had been 
delineated by hand calculations performed by Professor Lewis, This last 
reason was very important because it enabled the avthors to check their 


work and to ascertain the accuracy of the program, 
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The second problem programmed was the determination of stresses 
in a transverse web frame. A procedure similar to the one used in the 
Bureau of Ships was adapted to the computer. The numerous computations 
encountered in this problem made a computer solution very desirable. The 
calculations were made for a new destroyer type for which a partial check 
of the numerical results was available. Another important reason this 
problem was selected was to further the authors’ knowledge in structural 
analysis of a ship. The course of instruction in naval architecture had 


not included the calculation of transverse strength. 


PROPELLER SHAFT WHIRLING FREQUENCIES 


Description of the Problem 


Basically, the problem is to calculate the whirling or lateral 
frequency of a propeller shaft system in which the excitation is induced 
by the propeller. The method of calculation involves a numerical solution 


of a fourth degree differential equation of the form 


of which the general solution is of the form 


y = © sin kx + C, cos kx + &, sinh lx + C) cosh dx 


where 
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X=0 XK=L 


let y' = 8 and y" = ¥, Applying boundary conditions to determine the 


constants of integration: atx =0 
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From which 


s ow £OS 8 sinh B-sin B cosh 8B 


Lo sinh B-sin B 

s _. £98 cosh B-lL 

2 B( sinh B-sin 8B 

2B sin B sinh B 
°3 ° Sink Pesin B 
The above equations are applicable to one span of the shaft, 

In the solution of this problem these equations were used to determine the 
moment and angle of rotation at the end of one span and then carry these 
results over to the next span. The number of spans to be used is unspeci- 
fied, however, it is the opinion of the authors that four spans give 
sufficient accuracy. 


Since 


= 


various values of 8 are assumed and the residual moment at the free end 


e 
it 
ain 


plotted versus the frequency. The value at which the end moment becomes 
zero is a resonant point or natural whirling frequency of the shafting 
systeme 

In applying these equations it is useful to determine certain 
parameters in terms of the physical properties of the shaft system. These 


are given below and the notation is explained in Figure I. 





i bar 





Sign Convention 


Schematic Diagram of Shaft 
Figure I 


It can be shown from the equations previously developed that 
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polar moment of inertia of propeller 
and entrained water. 


the order of vibration being investigated, 
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The number subscripts refer to spans as shown in Figure I, 


From the above it can be seen that 


a 


P15 


where 


ark 


Here the second number subscripts refer to the span. 


be noted that 807 


sae 


lt 


- (34M, + C4453? 
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i should 


is assumed as unity. 


For determining Sj» Sys S4 for different spans it is convenient 


to use the following: 


where 
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The procedure described is repeated for each end moment and 


angle of rotation until the residual end moment, Ms), is determined. Then 


a new value of BS is assumed and the process repeated until a series of 


end moments versus frequencies results are generated, 


By plotting these 


values, or by interpolation, the critical whirling frequency can be found, 
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General Discussion of Program 


The method developed for solution of the propeller shaft 
whirling frequencies is ideally suited to computer solution because it 
is essentially one of numerical iteration, Various values of BS are 
assumed and the residual moment eaenated for any number of spans. %Since 
the labor involved in programming is the same regardless of the number 
of times a process is repeated, good accuracy can be obtained with no extra 
effort by the programmer and small increase in cost of operation of the 
computer, 

The underlying guide in prosramming this problem was to make 
it as general as possible. This was done so that anyone wishing to use 
this method in the future would merely have to change a few numbers to 
adapt this solution to his particular problem. The quantities to be 
changed were grouped separately at the beginning of the program and 
involve such factors as the physical characteristics of the propeller 
and shafting, the number of spans to be considered and the order of the 
vibration being investigated. 

The first step in programming was to assign floating addresses 
or "flads" to registers which were to hold the numerical quantities and 
key instructions. These flads contain the numbers or instructions upon 
which mathematical operations are to be performed. The purpose of 
floating addresses is to simplify programming by enabling the programmer 
to refer to a flad rather than an absolute address each time an instruc- 
tion is written, The advantage of this system is clearly demonstrated 
in the programs contained in the appendix, The flads used in pro= 


gramming and the numbers they represent are given in Appendix B, 
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In addition to storing the data and constants of any problem, 
registers must be assigned to store intermediate results and for temporary 
storage. 

The addresses of key instructions were assigned flads to aid in 
writing the program. For example, the flad, hl is the address of the first 
instruction in the program, After all the information has been fed into 
the computer the following instruction is used, "START AT h1" and the 
program beginS.e 

After the flad assignments had been made, the next step was to 
draw a flow diagram such as shown in Figure III, The purpose of this 
diagram was to aid in the organization of the program and to provide an orderly 
procedure for the programmer to follow, it was a picture of the problem 
showing the sequence of operations to be performed. The programmer then 
wrote a short routine that would do the calculations of each section shown 
in the flow chart. The last step, then, was to combine these separate 
routines into one program that would solve the entire problem. Of course 
considerable care must be exercised in writing and combining these routines 
Since the results of one section will influence the calculations in the 
following sections, 

When the values of frequency and residual end moment had been 
generated in the arithmetic element, they were transferred to a magnetic 
tape via an output procedure, When the program had been completed these 
results were then transferred to a flexowriter and typed. The results were 
presented in two columns of frequency vs. end moment. Thus, it was possible 
to determine the frequency at which the moment becomes zero by either 


linear interpolation or a plot of the data, 
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STRESS CALCULATIONS FOR A NON-SYMMETRICALLY LOADED WEB FRAME 
Description of the Problem 


A transverse web frame is a part of a system of frames and 
bulkheads that maintains the shape of the hull. A sketch of a web 
frame is given in Figure IVa. To do this it mst provide support for 
longitudinal strength members and carry vertical shear loads from the 
decks and hull bottom into the sides which are the main vertical shear 
carrying members of the structure. The web frame is a continuous beam 
ring that encircles the hull girth in a transverse plane, All loads 
applied to the frame are assumed to lie within this plane, but may have 
any orientation within that plane, Under such loading a web frame must 
resist axial loads, shear loads and bending loads. Since web frames are 
generally very long with respect to depth, bending moments are by far the 
most important aspect to be considered in the solution, 

Due to the complex interaction between the webs, longitudinals, 
and bulkheads, the exact solution of the probdlem of web frames and the 
system of members that maintain tne hull shape has defied rational engin- 
eering theory up to this time. Time limitations and the intent of this 
thesis prohibited an investigation of possible improvements in the method 

of the transverse strength calculations, It was decided to adapt to a 
computer solution a proven method, Although transverse streneth calcu- 
lations are not always performed for modern ships, the analysis normally 


(2) (6) 


utilized, if performed, is a long, laborious method developed by J. Bruhn 
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Figure IVb Web Frame Showing Redundant Forces 
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Bruhn's method is an adaptation of Castigliano's theorem of 
least work which states that "The redundant reaction components of a 
statically indeterminant structure are such as to make the internal 


(10) Determination of each of the redundant components 


work a minimum," 
requires that the differential of the total internal energy with respect 
to that redundant be equal to zero which is a minimum energy condition. 
Bruhn neglected the axial and shear components of energy and considered 


only the larger bending effects. Bruhn's basic energy equations for a 


web frame with two stanchions, as shown in Figure IVb, then become 
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Castigliano's first theorem states that, "the deflection in a 
certain direction at a certain point of a statically determinate structure 
is equal to the partial derivative of the total external work or the total 
internal energy with respect to a load applied at this point in the 
Fee Wcmorannendeiectianea<°)< Commereion iene etoticnliasindstecmmete 
structure into a statically determinate structure by means of cutting 
the structure and applying redundant forces as loads allows the application 
of Casticliano's first theorem, The above equations then represent the 
deflection between the limits of the integration. Since the web frame is 
continuous, the deflections and rotations at the point of action of the 
redundants Ms Po and dg are equal to zero, For the redundant colum 
forces R and S, the columns are considered imcompressible and the vertical 
deflection between the column ends is equal to zero. This idea of 
deflections and rotations is used in the description and analysis of the 
problem in Appendix B, 

To determine the redundants, five simultaneous equations must 
be formed and solved. Ina problem that contains more or fewer redundant 
supports such as stanchions or decks, the number of redundant equations 
required to solve the web frame would of necessity be changed to form 
one equation for each redundant force. 

After the equations had been developed and formed into a matrix, 
Crout's method of matrix reductionwis used to solve for the value of the 
redundants. These redundants were then combined with the applied loads 
and the moments recomputed. Utilizing the computed moments, the stresses 


due to bending were calculated, The direct forces were converted into axial 
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and shear stresses. These stresses and the bending stresses were then 
added algebraically to give the total stresses in the flange of the web 


frame and shell of the ship. 


General Discussion of Program 


The flow diagram contained in Figure V describes the sequence 
of calculations used in calculating the frame stresses, A cetailed 
description of the instructions utilized in each sub-routine and the 
sequence control program are contained in Appendix B. 

Due to the normal manner of construction of ships, the web frames 
are constructed as a symmetrical ring frame of variable cross-sections. 

The loads, however, are very complex and for a general solution must be 
considered as non-symmetrical, Note that for this problem only the simple 
static loads are considered. 

The contours of the ship's hull must be represented within the 
computer and, in order to find the neutral axis and to form the commonents 
of the loads, the slope of the contour must be computed at each station, 
The most obvious solution for the novice programmer is to form the equation 
of the contour and differentiate that equation to determine the slope. As 
novices this solution was attempted, but was soon abandoned due to the 
complexity of the program necessary to form the equation and the inherent 
error involved in first using approximate methods to form the equation 
and more approximate methods required to utilize the equation. The method 
finally selected was to choose a series of stations around the girth of 


the web and assume the web is represented by a series of short straight 
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Figure Vb Comparison and Alteration Routine 
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beam segments. Stations were selected in such a manner that the angle 

change between segments was both small and nearly equal which, due to 

the variable curvature, makes each segment unequal in length, A station 

was necessary at the point of action of each redundant. Such a proced- 

ure follows a common philosophy of digital computers in that where 

complexity must be weighed against a greater number of calculations 

and simplicity, the choice must be made in favor of simplicity. The 

number cf stations required to properly approximate the DD hull was 

h9 which should be sufficient to approximate the frame of any vessel, 
Sin @ and cos @ were determined directly from the coordinates 


by using the Pythagorean Theorem and basic trigonometric relationships. 


y va 2 : Ay Ax 
As =|/dx + day Sin @ = he cos 8 = As’ 


This was done for two reasons, first, only the sin and cos are used in 


Av 


the solution and secondly, the tan @ = ee 


becomes infinity when Ax = 0 
which is the case on ships with wall sides or tumble home. The comovter 
will not handle division by zero and the operations will be stopped when 
such a cas. occurs. To approximate the functions of the tangent angle at 
each station, the functions of the chord angles were averaged, While 

sin 6, + sin @,,/2 does not equal sin(@, + O72) the error will be small 
as long as the change in @ is small which is another example of using a 
Simple idea many times to obtain accuracy and generality. 


Using the total water head which in the present problem was 


the height of the main deck and the average of the y coordinates of the 
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ends of the segments, the water pressure loads were determined and vertical 
and horizontal load components formed. 
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Using the inputs of the web dimensions, the distance of the 
neutral axis (di, in ) from the hull shell was calevlated and the moment of 
inertia about the neutral axis determined — A. 


dike 


d ee 
NA, ‘ 


Ta. ™ teher1 ~ *eSu.a. 


Knowing dy , and the (sin ¢) and (cos g) the offsets of the 


neutral axis from the shell could be determined, 
xy a. = x~dy a (sin pd) 


YHA. = Ye-dy_, (cos hb). 


Using the same procedure developed for the shell, the segment lengths 


and functions of the angle were formed. 





27-6 


One difficult engineering decision was the consideration of the 
sharp corner at the junction of the side shell and the deck edge. In 
normal practice the corner is bracketed and represents a point of complex 
discontinuity which nearly defies analytical analysis, After a nurber 
of methods were considered, it was decided to make the corner a "hard 
spot" by inclucing the full bracket depth, The station spacing near tne 
corner was reduced to minimize the effect of the corner. In all cases 
the tansent at the corner was taken as -] which is equivalent to a slope 
angle of 135 degrees on the starboard side. This was done to form a 
smooth transition of the neutral axis around the corner, 

The moment cue to applied loads and water loads was deterrrined 
by ie H dy + Vdx. The moments due to the redundant forces were: 

Po + Gx + R(x-x, ) + S(x-x, ). 

The integrals of the moments were formed by the trapezoidal 
rule, While the trapezoidal rule implies linear segments between moment 
ordinates, this method was necessary to accomodate the random station 
spacing, tlhe error resulting from the trapezoidal rvle will be smell as 
long as the station spacinyz is close and the moment curve is smooth. Both 
conditions are present in this problem. 

After solving for the redundants anc recomnuting the moments and 
axial loads to include all forces, the cirect stresses and bending stresses 
were combined. The program as written cives the alternative of performing 
a check calculation on the set of frame scantlires or of alter‘ng the frame 
scantlings to arproach a more efficient desien, The stresses were compared 


with allowable stresses and the section altered by chancing the web flange 
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Figure VI Arrangement of Matrix Formed by the Integration Sub-routine 
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width by predetermined incremental amounts. If Qe aia OL) O, then the 


flange was reduced. If (ee ea < 0 then the flenge was increased. 

The outputs of each cycle change may be any commuted value the 
programmer desires, but in the present problem the outputs were flange 
width, stress at the flange and at the shell, direct stress and shear 
stress at each point. Other methods of scantling change could have been 
utilized or the scantlings could have been altered by a complete change 
of scantlings around the girth. In the latter case each cycle would have 
required a decision by the programmer, a procedure that has considerable 
merit where design experience in the mind of the programmer can be better 
utilized. In the present problem it was decided to use the conditional 
response as another example of the digital computer's capabilities. 

In the web frame problem the quantity of generated numbers 
necessary to produce a solution was so large that the auxiliary drum storage 
had to be utilized. Blocks of numbers were moved from the core memory onto 
the auxiliary drum to make room for calculations. Further, the size of 
the program was so large that optimum data output procedures could not be 
utilized without unwarranted complication of the program, The output data 
could have been typed or presented on films in columns of easily discernable 
numbers. In the present case the size of the program required output via 
an "fp" tape which gives five columns of numbers running consecutively 
across the lines and the address of the registers containing the numbers 
listed in the first colum, Judiciously piace zeros provide separation 


between numbers of different meaning. 
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To make the program as general as possible sub-routines were 
utilized to perform the calculations. A full description of the use of a 
sub-routine is contained in Appendix A. The sub-routines used will perform 
the calculation of stresses for a web frame of any dimension and shape 
desired. However, the number of redundants of the frame cannot exceed 
five and the water pressure head must equal the height of the deck edge, 
Due to limitations of storage space, the sequence instructions used will 
permit the use of only 9 stations around the girth, All of the above 
limitations can be exceeded by minor changes in the program should the 


necessity arise. 
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RESULTS 


The numerical results obtained from the computer olution of 
problem one are contained in Table I. Numerical results were not obtained 
for problem two. 

More important than the numerical results are the programs 
themselves. Problem one has been run and is known to be correct. The 
sub-routines of problem two, with the exception of the sub-routine for 
forming the stresses, were run and found to check, When all of the sub- 
routines were run in conjunction with the sequence control program, several 
minor programming errors were found. 

The first program will determine the critical whirling frequency 
of any propeller shaft where no more than four spans need be considered. 
It can be adapted to another shaft merely by changing the physical con- 
stants of the shaft involved, 

The second progran will, if the above mentioned errors are 
removed, determine bending, axial and shear stresses in a transverse 
web frame and, if desired, will adjust the size of the web frame until 
the computed stresses approach the allowable stresses. The program is 
general within the limitations stated below; the required inputs being 
the dimensions of the web frame, the location of the stanchions and 
coordinates or points on the transverse section, the weight and shear 
loads at these points and the hydrostatic design head above the keel, 


The following are the limitations of the program: 





ay 
(2) 


(3) 


(1) 


) 
(6) 


the frame scantlings are symmetrical about the ship centerline; 


the points defining the transverse section must be symmetrical 
about the ship centerline and only 25 may be taken on one 
sides 


one of these points must be on the centerline keel, one 
on the deck centerline and one at each end of all the 
stanchions $3 

the hydrostatic pressure head must be equal to or greater 
than the height of the deck edge and hydrostatic pressure 
loads are not generated for the decks 

ship must be in upright position; 

only one centerline or two off-center vertical stanchions 


may be used. 


Computer time for calculation of problem one was 10 minutes, 


but the typing time of the results was 20 minutes. Computer time for 


the solution of problem two was 10 minutes, 





Frequencies from 
hand calculations 


+6 ORDER 
3.69 cps 


ee ial 
8.20 
11.9 
15. 00 
21.00 


~6 ORDER 


Not Available 


+1 ORDER 


Not Available 


336 


Frequencies from 
computer solutions 


3.870 cps 
5.656 
Beoes 
11,762 
1 Seie39 
200085 


B00 
5.656 
2 ea 
12702 
329 
208915 


3.870 
5.656 
8.527 
Tae (Oe 
15.239 
20.915 


Critical Whirling Frequencies for Propeller Shaft 


Table I 
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During the period that the thesis was being printed, the pro- 


gramming errors in the second program were remedied and the results listed 


in Table Il were obtained. 


Table II contains the values of the redundants 


end the stresses at the keel and the deck at centerline, Due to the 


difference in station spacing between hand and machine calculations, the 


centerline points are the only coincident stations for the two methods 


or Calculation. 


R 
S) 
om = p/A 


Trlance keel 


J shell 

Ta 

Jp eee 
vs 


Hand Method 
47/250 aie. 
+20, 81 
(assumed zero) 
41925 


(assumed equal to R) 


~3880 */sq.in. 
+30,750 #/sq.in. 
“15,40 */sq.in 
+1.70 */sq.in.e 
-10, 80 #/sq.in. 


+3, 350 #/sq.in. 


Table II 


Computer Method 
416,395 -fte 
+19, 805 

+282 

+662 


+662 


-5193 #/so.in. 
$10,070 “/eq. ii. 
-16, 586 */saeine 
C00 i Saeane 
+25,000 */sq.in 


-12,868 */sq.in. 


Results of Web Frame Stress Calculations 
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As shown in Table II there is little correlation between 
hand and machine calculated values, However, this lack of comparison is 
caused by two conditions. The most important condition is that the 
hand calculations have been found to be of doubtful validity. ‘Several hand 
calculation errors were found and the scantlings used are not defined, 
Recent information indicates that the scantling valves used in the 
Bureau of Ships hand calculations do not correspond with the scantling 
values presented in the contract plans used as a basis for machine 
calculations. 

The other condition that causes a small error is the fact 
that summation of vertical forces does not equal zero. The difference 
between weight forces and buoyant forces must be balanced by vertical 
Shear forces applied along the vertical sides of the hull. This shear 
force was taken from the hand calculations and was included in the 
input values for the program, The machine computed buoyant forces were 
different from hand computed values due largely to the difference in 
station spacing and the method of averaging the water pressure. Therefore, 
the input values of shear did not balance vertical forces computed by 
the machine. It is now obvious that the shear force values must be 
computed by a sub-routine in the computer program. 

This lack of the means of proper comparison points out one 
other general rule of programming. The programmer should perform hand 
calculations for the problem being programmed using the method adapted 


to the computer. The authors decided that the use of existing hand 


~ +? 
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computations would be a time saving device; however, as the results 
indicate, this was a poor decision and it is hoped that future prorrammers 
will benefit by the author's experience. 

With the exception of the shear error discussed above, the 
program does work and produces good results. With slight correction the 
web frame program is suited for research and ensineering calculations. 

The time required for machine calculations was 25 seconds; 
however, program conversion required five minutes and typing time for 
results of calculations required twenty minutes. 

Time limitations did not permit the running of the comparison 


and alteration program for the web frame, 
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DISCUSSION OF RESULTS 


The subject of computer application to naval architectural 
problems is as general as the field of naval architecture itself. Each 
proposed application will require a separate analysis of feasibility. 
However, certain general conclusions and impressions can be gained from 
the two problem examples of this thesis. To each general conclusion there 
exists certain possible exceptions; however, these exceptions only tend 
to prove the contention that digital computers are not universal panaceas, 
but only local aids, 

The programming errors found in problem two are normal errors 
that occur in every program and do not reflect upon the method of approach. 
The nature of these errors are incorrect addresses, typographical mistakes 
and human errors. The correction of these difficulties is only a matter 
of time and patience. 

The results of the two progrems indicate that digital computers 
are applicable to long and complicated naval architecture problems and 
that the numerical results are more than comparable with those of hand 
calculations, Of greater engineering importance than accuracy, was the 
great reduction jn calculation time obtained by computer usage. For the 
first problem the computer required 10 minutes as opposed to an estimated 
hand calculation time of thirty (30) hours and in the second problem a 
more defined superiority is estimated in that the computer solution will 
require approximately ten (10) minutes versus an estimated sixty (60) 


hours by hand. However, this comparison does not depict the true time 
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picture and the long tedious hours of program writing and ever present 
"debugging" or program correction must be considered. Although the pro- 
gramming time varies greatly with the programmer's familiarity with the 
problem, computer and also, with the nature of the problem itself; in 
general the programming time and the time required for one hand calculation 
are roughly commensurate. The authors estimate that a practicing naval 
architect very familiar with the proposed problem and experienced in 
programming could program problem one in twenty (20) hours and problem 
two in forty (10) hours. 

The "debugging" time for each problem would require another 
ten (10) hours. The fact that the authors, unfamiliar with both the 
problem and computer, required six times that amount of time indicates that 
experienced personnel are a very important application factor. In the 
present case problem understanding and organization were responsible for 
the greatest loss of time. As a general statement it can be said that for 
an engineer to properly program a particular problem, his knovledge of 
that problem must be more thorouch than that required for hand calculations, 
Not only must he know the individual steps in the problem, but also must 
predict the effect of each step upon the following step and the general 
value of the numbers generated. 

Also, the authors found that one programmer workins alone is 
just as effective as several working tovether on the same problem. The 
reason for this is the necessity for a program to follow a single trend 


of thought and the exact accounting required. 
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To provide a more comclete time description of pvrogramins, 
the nature of the proposed problem should be analyzed. Problem solutions 
that involve repetition of the same equations will tend to have a much 
shorter programming time. For instance the time required to program 
problem one for lO spans of shafting would be only slightly longer than 
that required to program spans. Also, the program required to desig 
the web frame by cycling and altering the scantlings is little longer 
than the program for check calculations. 

After the solution of one particular problem by a computer a 
general program will find further use in solving other similar problems. 
The programs written for this thesis are reasonably general and the 
limitations stated in the results are not of serious damage to the program. 
The limitation of problem one to four shaft spans in a reasonable engin- 
eering division in that propeller excited whirling vibration amplitudes 
beyond span four are so small that their effect on the end fixity of 
span four is insignificant. In problem two the limitations were imposed 
to save storage space and reduce the complexity of the program. ‘This 
coincides with the authors! opinions that for engineering applications 
several simple programs are easier and faster to handle and program than 
one complex all purpose program. For instance, web frames are normally 
constructed symmetrically, 25 points should properly define most normal 
web shapes and to have a station at each major load or redundant is 
reasonable. The design hydrostatic head is normally above the maindeck 
and waterloads are normally considered a constant value along the deck, 


While the listed position of the ship is an important consideration, 
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modern naval architects consider the upright position for practical 
calculations and should the listed condition become necessary for re- 
search calculations, rotation formulas and non-symmetrical water loadings 
must be incorporated into the program. To change stanchions and add 
decks as redundant supports major alteration of integration procedures 
would be required. 

A general program is useful in design and research to compare 
and generate great quantities of data and to test the results of slight 
alterations upon a major problem. For instance, program one covld analyze 
and compare the whirling frequencies of hundreds of ships in a relatively 
short time or the effects of propeller weight, fixity and blade number 
may be determined for numerous combinations of values. Problem two 
would find use in generating design schantlings for various load and size 
variations such as those found in classification society publications, 

As a sub-routine in conjunction with a more complex program the stiffness 
of web frames relative to transverse bulkheads could be determined by 
application of the theory of beams on elastic foundations. By reiteration 
of calculations under varying loads the variant or impective loadings 

may be analyzed. Since the longest time consumed in a program is in the 
program writing, problem repetition increases the efficiency of the com- 
puter application. 

To be feasible for engineering calculations the computer must 
save time, money and increase the effectiveness of the engineers it aids. 
Therefore, the computer must accent numbers of enrineering magnitude and 


quantity in a convenient form. The coded instructions must be straichtforward 
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and have names descriptive of their use. Also, the generated results must 
be presented in such a manner that wide spread publication is possible. 
Whirlwind I and some commercially available machines will accept decimal 
digits of sufficient length and quantity to perform most practical 
engineering and research problems, Of extreme importance is the necessity 
of maintaining a large library of sub-routines capable of performing such 
calculations as sin ~@, \/X, etc. Frequently it is desirable to purchese 
the services of an experienced programmer to act as a consultant when 
first using a computer. 

In computer application there exists two basic schools of 
thought. One group thinks that problems should be programmed by experts 
who perform no other function than prorramming. In this case the engineer 
presents the problem to the programmer and at some later time is given a 
solution. The other group states that a computer is best utilized by 
having engineers familiar with the problemeffect the programming. While 
there exists a middle road somewhere between the two extremes, the avio™ — 
are of the opinion that engineers should program the computer. This 
performs two useful functions for the profession as a whole. One, the 
programmer can make the necessary engineering decisions to adapt the 
problem more successfully and secondly, the engineer becomes more aware 
of possible computer epplications in his field. 

The qualities required of a programmer are no more exacting 
than those required of a food naval architect. He mst have almost 
infinite patience, imarination, ingenuity, sood ability to organize 


problems and a sood background in both theoretical and numerical methods 


eo 
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of mathematics. The work is tedious and all instructions must te perfect 
Since the machine will scenerate wrong answers as fast as correct ones, 
Digital computers are often referred to as "brains," but nothings could be 
more incorrect. in that thought is impossible. Actually a computer more 
nearly resembles a desk Calculator or an abacus which operates at a hich 
rate of speed. Imagination and ingenuity are so necessary to produce a 
workable, general program that many commercial computing companies obtain 
their programmers from chess clubs, cross-word puzzle enthusiasts and 
other hobbies where deep concentration and original thought are required. 

The future of digital computers in naval architecture appears 
very promising from both the research and practical viewpoints, In 
research they can be used for the rapid assimilation and processing of 
large quantities of technical data and for the solution of problems that 
require an inordinate amount of time to calculate by hand. Tvpical examples 
of the d ata processing type of problem would be a comparison of stability 
among several different vessels and other similar problems involving 
statistical analysis to establish design criteria. Since the same pro~ 
gram may be used to solve a number of similar problems, design data in 
the form of charts and tables can be generated for the grillare beam 
analysis of decks and other complex problems. such a procedure permits 
the determination and testing of important parameters for later anplica- 
tion to specific designs, 

For specific design applications comouters may be used to effect 
large savings in time and labor, This is possible because the computer 


can approach more closely an optimum design with little increase in 
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computing time and expenditure of effort. In this category of aovli- 
cation are midship section design, transverse framing design and flight 
deck bents. In addition, refinements intheory end in the extent of calcu- 
lations can also be utilized to a greater degree than practical by hand. 
Examples of these refinements are the use of prillage beam analysis for 
decks and bulkheads, the performance of stress calculations for a ship 
in drydock, design of propellers by circulation theory, and the investi- 
gation of submarine dynamical stability. In addition to desien calcu- 
lations a digital computer can be applied very readily to scheduling 

and accounting problems. In shipbuilding this system can be used in 
determining the weight and center of gravity of a ship as it is being 
built. Another possibility is merchant ship carro arrangement, 

Theses present an opportunity for wide and varied use of 
digital computers. In this application the comouter should be need 2s 
an aid and not as the general subject. A student attempting a problem 
in which the magnitude of the calculations is the dominant factor should 
consider computer aid. A desireable arrangement is to have a team with 
one member learning the problem and another member learning the comouter, 
In this manner both can derive maximum benefit from the thesis and a 
larger contribution to the profession can be made. Since computers are 
becoming more widely used, familiarity with their capabilities may be 


of considerable future value to the student, 





CONCLUSIONS 


Digital computers are applicable to naval architecture problems, 


The average naval architect has sufficient background to prorram 
a digital computer. 


Some general requirements of a programmer are: 


ae They should have a complete engineering understanding 
of the problem to be programmed. 


b. Familiarity with numerical methods of mathematical 
processes is necessary. 


c. They should be amenable to tedious work. 
The proolem to be programmed must be suitable to comcuters since 
it may take longer to program than to determine one solution by 
hand e 


One programmer can work faster than two or more. 


In general the limitations imposed on a problemare not serious and. 
are no more restrictive than those in hand calculations, 


Complete generality of a program is not always feasible. 


For engineering applications the following features of computers 
are necessary: 


a. Computer word-length should be able to handle engineering 
numbers. 


b. Library of sub-routines should be available. 


c. Computer should be easy to program. 





RECOMMENDATIONS 

Digital computers should be used to solve naval architectural 
problems. 
Digital computers should be used in theses and research problems. 
Any future thesis should consider the computer as an aid in the 
solution of a specific problem and not applicability of computers 
to a general field. 
For engineering applications the programmer should have; 

a. good engineering understanding of his field; 

b. knowledge of numerical methods; 

Cc. organizational qualities and a great deal of inseenuity; 


de. aworking knowledge of the computer utilized. 
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Supplementary Introduction 


A computer program is a series of logical instructions 
telling the computer what mathematical operations are to be performed 
on the data. Thus there are two proups of inputs, the data and the 
instructions. The data fed into the computer is stored in the core 
memory. Hach number or piece of data is stored in a register or "cell" 
the location of which is designated by a number called the address of 
that register. After all the data is stored the instructions are 
stored in the unused portion of the core memory. The instructions are 
stored in sequence and will be performed in the order in which they 
are stored. No unused registers are permitted to interrupt the storage 
as these will stop the computer. The programmer may wish to know the 
address of key instructions such as the one that starts the computation 
process, He may assign an address to these instructions or allow the 
computer to do so, as explained later, 

The instructions given to the computer consist of two parts, 
the operation and address sections. The operation section is first and 
uses letters to tell the computer what mathematical operation is to be 
performed. The address section gives the address of the register contain- 
ing the data on which the operation is to be performed. The instructions 
recognized by Whirlwind I are listed in the glossary. 

When writing a program for a digital computer, there are several 
useful devices that simplify the procedure, The programming devices described 
below apply to Whirlwind I and may or may not apply to other computers. 

Probably the most useful device is a floating address. A 


floating address or flad is used to eliminate the necessity of referring 





to the absolute register address in the core memory when writing the 
program. At the beginning of a program the inputs or key instructions 
on which mathematical operations are to be performed are given floating 
addresses. These addresses consist of a letter and a number, for 
example al,a2 etc. Any letter except O and L may be used and the sum 
of all the numbers behind the letters may not exceed 255, In writing 
the program, then, the flad is used in the address section of the 
instruction to represent the absolute address. The computer will assign 
absolute addresses or register numbers, to these flads unless this is 
done by the programmer. For example, if al is the first flad in the prog-~ 
ram, the computer will assign it to register 32. If, however, the prog~ 
rammer wants al in the register 130, he writes 130/al and the computer 
will place al in register 130. A table is usually produced by the 
computer showing absolute addresses assigned. It should be remembered 
that a flad can have only one value in any one program, 

Another useful device is preset parameters. These are numbers 
other than data that are used in the program. Preset parameters are 
frequently used to designate the number of times a computation is to be 
repeated. They differ from flads in that they may have any absolute 
value desired whereas flads must have the absolute address of a register. 
Preset parameters must be assigned a value before they are used and they 
Will retain this value until it is specifically changed. Preset parameters 
are designated by two letters and a number. The first letter must be 
P, u or 2 and the second letter and the number may be anything, except the 


letters o and 1 cannot be used, 
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In writing a program for a digital computer subroutines are 
very helpful. Basically there are two distinct types of subroutines which, 
for lack of a better name, shall be called internal and external. /JInternal 
subroutines are those performed automatically by the computer upon receipt 
of a single instruction. An example of an internal subroutine is the 
group of output instructions listed in the glossary. Each of these 
instructions will cause the computer to go through a certain predeter- 
mined procedure without any further information from the programmer. 
External subroutines are further divided into two categories; one category 
is designated library subroutines and are perfectly general for any prog- 
ram. The other external subroutines are groups of instructions that 
perform a repetitive section of a particular program, Library subroutines 
are groups of instructions usually written by an expert programmer for 
the purpose of performing a particular mathematical operation, They 
may be an integral part of the main program or they may merely serve as 
an auxiliary program. Examples of these external library subroutines 
are groups of instructions used to determine the square root of a number, 
trigonometric and hyperbolic functions of angles and the solutions to 
Simultaneous equations, Most computer installations have libraries of 
subroutines available for programmers. These are written quite arbitrar- 
ily by using relative addresses and can be used with any main program, 
There are some rules that must be followed when using subroutines. 
The main program must be interrupted to call the subroutine into use. 
After the subroutine has performed its function, it mst direct the con- 


trol element back to the proper instruction in the main program, An example 
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of a subroutine follows: 


Major program 
ia ae 
£6, isp g3 -- This instruction directs control to follow the 


instruction in the register designated by g3. 


SUB-ROUT INE 

¢3, ita g3 plus 20 ~-- This instruction 
places the address contained in the 
control element in g3 plus 20 and adds 
one to that address. In this example 
the address contained in the control 
element is f6. 

Tadeo 


more instructions 


g3 plus 20, isp f6 plus 1 -- The address of 
the isp instruction was supplied by the 
ita instruction at the beginning of the 
subroutines 
16 plus Iy iad a3 
The authors do not wish to imply that a program should consist 
of a series of subroutines, but they are very useful in programs that 


perform similar operations many times throughout the problem, 
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The numbers usually encountered in engineering are too large 
to be stored in one register which, it will be recalled, is 16-binary 
or about )} decimal digits long. To accommodate these larger numbers 
a technique called floating point arithmatic is used. Two consecutive 
registers, giving 32 binary digits available for each number, are used 
to replace the single register. Some of these 32 digits are used to 
record the actual number and some to locate the decimal point. In this 
manner the range of significant digits can be increased by restricting 
the number of digits used to locate the decimal point. In the binary 
system the decimal point is written as some power of 2 and usually 6 
binary digits are used to express this power. Since the first digit 
of each register is used for the sign of the quantity in that register, 
in this case the sign of the second register is the sign of the exponent 
of 2. Twenty-four digits are left to express the significant number. This 
particular distribution of binary digits is referred to as the (2.6) sys- 
tem and can handle numbers equivalent to eight decimal digits with expon- 
ents up to 1073, Other systems are available and the size of the numbers 
in the problem determines which one should be used. 

In order to manipulate these two register numbers, an inter- 
pretive mode has been internally programmed to function as a part of 
Whirlwind I. All instructions that refer to double length registers 
are preceded by an i. If calculations or operations are to be performed 
on single length registers, the same instructions may be used without the 
i preceding. Hence there are two modes of calculation; interpretive 


for double register length numbers and whirlwind for single length register 
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numbers. It is possible to pass from the whirlwind mode to the inter- 
pretive mode by use of the instruction IN and to reverse the direction 


by use of the instruction OUT. 
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GLOSSARY OF INSTRUCTIONS 


TERM DEFINITION 


eax clear MRA and add contents of 
register xX. 


ics X clear MRA and subtract contents 
of register x. 


iad x add contents of register x to 
what is in MRA. 


isu xX subtract contents of register x 
from what is in MRA. 


imr x multiply what is in MRA by the 
contents of register x and round 
off product to fifteen digits. 


aie higae.< divide contents of MRA by the 
contents of register x. 


its x transfer the contents of the MRA 
to register x losing what was 
previously in register x. The 
MRA remains unchanged. 


iex x exchange the contents of the MRA 
with the contents of register x. 


isp x take the next instruction from 
register x. Does not affect the 
contents of the MRA® An uncon- 
ditional instruction used to break 
the sequence of operations. 


Tep x if the contents of the MRA are 
negative, take the next instruc- 
tion from register x and continue 
from there. If the contents of 
the MRA are positive, take the 
next instruction in seqence. A 
conditional instruction used to 
break the sequence of operations. 





Glossary of Instructions (Continued) 


TERM 


oc. 


Ver xX 


ict x 


ici x 
ied x 


iti x 


Ja xX 


ita 
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DEFINITION 


select counter number x. Without 
this instruction counter zero 
will be used whenever +C appears. 


cycle reset. Sets index register 
of counter to zero and criterion 
register to x. 


cycle transfer. Increases con- 

tents of index register by one. 

If contents of index register 

greater or equal to contents of 
criterion register, set index 
register to zero and do next 
instruction in sequence. If con- 
tents of index register less than 
contents of criterion register, 

take next instruction from register x. 


increase the contents of index 
register by number x. 


decrease the contents of index 
register by numer x. 


transfer the right eleven digits 
of the index register into the 
right eleven digits of register x. 


add contents of index register to 
the contents of register x and 
store the result in the index 
register and register x. 


replace the address section of the 
instruction in register x with the 
address that is one more than the 
address of the register containing 
the last isp or icp if the contents 
of the MRA are negative. 
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Glossary of Instructions (Continued) 


TERM DEF INUTION 


i TOA record the contents of MRA on 
direct printer (typewriter). 


{MOA record the contents of MRA on 
magnetic tape for delayed printing. 


150A record the contents of MRA on 
oscilloscope for photographing. 


iFOR this instruction provides an 
automatic device for obtaining a 
suitable layout of output data 
in columns, lines or blocks, 


Note: The above output instructions are usually followed by a series 
of letters and numbers that indicate the desired form and 
arrangement of the output, 





Appendix B 


Details of Procedure 








526 


Details of Program for Determining the Whirling 
Frequencies of Propeller Shafts 


The first step in programming this problem is to put the 
numerical values of various quantities into the computer. This is done 
by the use of floating addresses, or "flads."" The proper use and func- 
tion of flads is explained in Appendix A. 

The following table lists the flads, the quantity and the 


values used. Figure ]| defines the symbols and nomenclature. 


ail L, = 77" - length of tail section. 
a2, L, = 513.6" ~ length of inboard sections. 
L, = 576.0" 
‘ee ou. 6" 
3 
Di. Bato = 17el} ~ 10°" on - E times moment of inertia of 
. tail section. 
D2; el =) 85} We ome 1b-in® - E times moment of inertia of 
inboard sections. 
BI, = 12.4 x 10° Ib-in® 
EI, = 12.4 x 107° ip-in® 
cr, i = 17,815 lb ~ weight of propeller 
ore Vy = 28h) If in? - specific weight of material used in 


inboard shaft sections, 
Yo = «28h 1bf in? 
Yq = +28 Lof in? 


adi; AS = 19325 in’ - average area of shaft cross section 
for each span. 


A = 17.5 in’ 


A = eS iné 
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az. 3 - number of shaft spans used. 


as; 135 - number of increments of B desired to 
be used. 8 is increased by increments 
of 2 degrees. 


6 


y 


1 
eo, 3 sa q 
q is the order of frequency being 
25) _ 6 — investigated. Plus sign is used 
q pee for backward precession and the minus 
sign for forward precession, 


rae q=il 
eh, J = 2hh.0 % 10° - polar moment of inertia of 
propeller and entrained water, 
eb, BS = 90 - initial starting value of 8,5 in degrees, 


The above listed flads are the only ones that must be changed 
in order to adapt this program to another problem, In order to accomplish 
this change it is merely necessary to add a tape with the correct values 
at the end of the program. 

The following table lists flads that contain constants which 


will remain unchanged for every problem: 





els To = 3,116 
e2, 180 
e5, g = 386 in-sec* 
e7. BS = current value of B, (8, + increment) 
yh 
a? 180? 
Tr 
eo. 
(180)* 
el0, 2 
ell, 1 
el2, 3 


el3, 4k 
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A block of registers consisting of f and g flads has been 
set aside as intermediate storage. This is necessary to enable the 
computer to calculate quantities, store them temporarily and use them 
later in the program, Another group of registers consisting of t and u 
flads are set aside for temporary storage, These quantities are used 
once and then discarded. 

The last group of flads used is the h and m group. These are 
used to control the program by assigning key instructions to them, 

The outputs from this program will consist of a series of 
frequencies and corresponding end moments, These values will be arranged 
in two columns with the frequency corresponding to the moment that 
appears on the same line. lt is very important to include the sign of 
the moment, for we are more interested in the point where the moment 
passes through zero than the actual values obtained, 

In performing certain routine arithmetical operations, such 
as finding the square root, use was made of the sub-routine library in 
the computer laboratory. These sub-routines are very reliable and save 
considerable time when programming. 

The actual orders used for programming this problem follow. 
The short paragraph at the beginning of each section explains briefly 


the function of the orders in that section, 
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These three orders place the sub-routines for square root, 
sine-cosine and sinh-cosh at the beginning of the program, This is 


very desirable as it eliminates any confusion with the main program. 


Order 
Sl, LSR, FU8 
S2, LSR, FUha 
S3, LSR, FU5 


Put in flad assignments al through glO as previously noted, 


These orders are assigning values of 3, 2, 3 in the flads 
Mi, M2, M3, respectively, the function of these flads is to set the 
cycle counter criterion registers to a predetermined value. The opera- 
tion is done in WW] since the numbers are single register numbers. 


The last order is used to enter the CS mode. 


Order Input to MRA in MRA Transfer 
from MRA 
ald Oa oe 8 3 > 
td, ml 3 
td, m3 S 
Sune ale 2 
td, m2 2 2 
IN 


This group of orders will calculate the weights of each shaft 
span and store the results in fl. The first order says that the opera- 
tion is to be repeated three times. The +0 at the end of the next 
four orders is to permit the address section of the instruction to be 


increased by two registers on each cycle thus allowing three different 
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weights to be calculated by the same set of orders. The last instruc- 
tion is necessary to return control to the first instruction at the end 
of each calculation, to increase the address of each instruction oy two 
registers after each operation, and to pass control to the next order 
when the operation has been performed the number of times specified by 


the first order, 


Order input, toehink In MRA Transfer 
from MRA 
los eg 
h2, ica, c2+c 15 Ve 
imr, a2+c x L, ae A 
imr, dl+c x A, L, x ¥, % A, 
its, 11+¢ W, W, 


This group of orders will calculate - and store the results 
in f2. The cycle will be performed three times since the order ml, icr 
is still effective. 


Order Input to MRA In MRA Transfer 
from MRA 





hy tea, bete EL BL, 


E_ I 
idv, a2+c Sa — 
1 L 
1 
its, f2+c a g 


ict, h3 





Bae 


This group of orders will calculate C, for each span and 


d 


store the results in f3. The cycle will be repeated three times since 


thevorder mi. acr 15 stiil effective, 


Order Input to MRA In MRA Transfer 
from MRA 
hiaica, t2+c — Se 
Bt, 
TGV fae +e — Ly 3 
EI 
idv, a2+c +L, — 
EI 
idv, fl+te ih = 
m1) 
1 
EI 
its, £3+c aoe C. 
W.L 
ae 
ict, hh 


This group of orders will calculate the constant /(180)° 
and store the result in e9, This operation will only be performed once 


as no cycle is indicated. 





Order 


ica, 
idv, 
idv, 


its, 


el 


e2 


e2 


Input to MRA 


— 180 


— 180 


Gos 


In MRA Transfer 
from MRA 
tT 
7/180 
2 
1/(180) 


n/(180)* /(180)° 


This group of orders will calculate the constant (1 /180)" and 


store the result in e8, 


The device used here is that when a number is 


transferred from the MRA to a register the MRA retains the number. The 


fourth power is formed by multiplying the number in the MRA by the number 


that was transferred which is the same as squaring the number in the 


MRA, in this case (1/180)*. 


Order 


imr, 
bss 
imr , 


its, 


el 


e8 


e8 


e8 


Input to MRA 


x (1/180)* 


In MRA Transfer 
from MRA 
(1 /180)° 


(n/180)*  — (#/180)° 
(n/180)4 
(1/180)4 (n/180)4 


This group of orders will calculate os and store the result 


in fh. 








Order Input to MRA In MRA Transfer 
from MRA 
eile let 
Za O72 ee 3 See 
ied rig 
WE 
jag heme | x W pia 
Pp W.L 
el: 
WLEI 
imr, al x L polit 
=! O W 3 
le 
W LE . 
imr, al x L L 
O W 3 
4 
imr, al x Ly 3 
Soak 
8 
W LAE TL 
idv, bl aT pee EL 
oe W 1B i 
1 00 
WL? 
imr, 8 co “pod (4): 
180 W 35 ir 180 
4 OO 
S 
W L-E 
nee) "POL yt Fe 


3 180 
rigaheoue 


This group of orders will calculate C. and store the result in 
£6. In this operation the square root sub-routine is used. The third 
order is used to enter this sub-routine, which is in the Whirlwind mode, 


requiring the IN order to return to the CS mode once the square root has 


been calculated, 





60. 








Order Input to MRA In MRA Transfer 
from MRA 
ica, e5 g g 
E_I Ig 
. ea a 
a (i % a 3 
Mil ariel 
E_i¢ 
SD, sl] ae 
L 
Wy 1 
a ale: 
L 
vr 1 
E 
: Wt , Hee i448 
eS X (180) (180) WID 
a 
idv, e10 3 a? rah 
4 ; 2(180) WL2 
aes 8 





| peer a 


This group of orders will calculate Fe and store the result 
in f7. The operation is performed twice, since m2 is 2, because this 
constant is used to determine B, and By. If more spans were used this 
operation would be performed a number of times equal to two less than 


the number of spans. The square root sub-routine is used twice in 


succession here to form the fourth root of a number, 





aie 


Order Input to MRA In MRA Transfer 
m2, 2cr 
E 
: ah 14 
Ho aicane £3 5) as 
Li 1 set 
BT, B Ty4,15 
idv, £342+c — eS) ae 
EL: bee 
Wo 2 Tl 5 2 
sp, sl 
sp, sl 


ioe eae 





ict, h5 


This group of orders will calculate ie and store the result 
in f5, This operation is to be performed three times as indicated by 
the first order, The flad ml was not used here because the number of 
times this operation is performed depends on the number of vibrational 
orders being investigated and is independent of the number of spans, 
Usually three orders of vibration will be investigated, (synchronous, 


forward and backward precession). The flad h6 is the beginning of 





the first loop of the program, This means that when the program has 
been run through once it will automatically come back to h6 and do the 
entire program again using the second value of e3. This loop will be 


performed three times as indicated by the first instruction. 








Order Input to MRA In MRA Transfer 
from MRA 
ab eh are. 
E, 
[ole Pisaike ts (aps a: =a a 
W 
ila age 
E_IJ 
imr, el Xx J 1 
ea 
EL JL 
imr, al XL, 1a 
W.L 
el 
T_JL 
idv, bl eal “4 = 
ET WL 
OO 
IJ 
imr, e3+c x (S44) L190 2, 4) 
q Twi-* 4 
OF Oil 
| eae Kier a 
. Sel! (site Oe alia i matted 
Og 6 Sal age 
K. I. JL 
; gy ee 0 ae ee eae 
ts, f Et HS) 
ts, 1? nae % * 4) FB? oe 
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This group of orders assigns a value to the flad mj. The value 
used is that contained in d3 which for this problem was 135. This number 
was selected because we wanted 8 to vary from 90° GO 360° in as small an 
increment as possible to increase the accuracy. Due to limitations on 
the magnetic tape output the smallest increment feasible was on The 
first instruction is used to select the Whirlwind mode since the trans- 
fer of digits as indicated by the third instruction cannot be performed 
in the cs mode. The fourth order is used to return to the cs mode 
before proceeding with the remainder of the programe 

The last instruction is used to select a new cycle counter, 
This instruction permits a loop within a loop to be performed. In the 
preceeding group of instructions counter zero was used and is stepped 
once each time the sequence is performed, but a different counter is 
used to allow the following sequence to be performed, and counted, a 


specified number of times without interfering with the previous counter. 


Order Input _to MRA In BRA Transfer 
from MRA 
OUT 
Cad 135 135 
td, m4 es 135 
IN 
Sy vet 


This group of orders will calculate the frequency and print 


the result on the magnetic tape output device, This operation will be 
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performed 135 times, indicated by the first order, for each cycle of 

the main program. The second through the fifth orders calculate Bo 

and git and store them in gl and g2 respectively. The frequency is then 
found by miltiplying OP by Bae The last instruction is an order to the 
computer selecting the output device and giving the information necessary 


to print the frequencies in a colwm, the number of digits desired, etc. 











Order Input to MRA In MRA Transfer 
from MRA 
mi, icr 
h7, ica, e6 90° 90° 
imr, e6 x 90° (90°)* 
: 0,2 
LUS, gl (90 ) B., 
imr, gl vaneoe (90°) 
its, g2 (90°)! By, 
6 mn Ey he ee ae 
LCay (18 ye WL? ; 2 12 
2(180 yy 2(180) Ww, Ly 


n(90)* [*148 


e 3 
2(180) WL, 





imr, gl % (90°)¢ 


imoap 123).56ssss 


This group of orders will calculate the moment at the inboard 


end of the tail shaft, Mm, 


will also be performed 135 times for each cycle of the main program. 


12 and store the results in g5, This operation 





Cos 


p. and Dy both functions of 8 are comouted first and stored in g3 and gl 
respectively. Using the forma in the flow-chart they are combined to 
calculate M4 It should be noted that this and the preceeding calcula- 
tions do not require more than two double lengtn registers to store the 
results. tn the frequency calculations the results are printed out as 

soon as they are obtained. Here the results stored in g5 are used immedi- 
ately in the subsequent calculations after which they are no longer desired. 


Therefore, the m., calculated is lost as each succeeding moment is determined. 


OL 


Order Input to MRA In MRA Transfer 
from MRA 
5) S 
W L°-E.I W L-E.I 
ice th 0 POLL Gy BOR Tl 
W.L-E.t. 20 Wirz 10 
i 00 He OO 
W LE 
hy nee mB \)y 
Wiigened: 
pF oo 
W LE 
Spiel 7B OY 
its g3 ) p 
Wier 180 w 
OO 


fan fo “Ao 4 by yk Aa” 2. yy 
erwi3?2 @ 10° pry 732 4 180 


OO: bo O-Onr 


BLL JL 
slam +L, 


imr e6 B ty kt 
SFr oleate 2) 
oobi 
E_LodJL 
its gh ao a +B By 


32 
“oo 





Order 


ica g3 


idv ele 


bos +) UL. 


18 hier s 1S, 


jad lh9 


isu tl 


isu gh 


VES ee 


Lea) “tL 


isu gh 








Input to MRA In MRA 
(continued) 
oe ae 
£3 fy 
5 
he 
5 
a ee 
slip 
AP 
* fg 2 
an 
+ 1 19 
,, APs 
3 cies 
one 
A, “fT? 
fps 
2 
ee p,, 
37 3 
- By p ys 





Transfer 
from MRA 
fy 
3 

an 
3 
fay 
3 J 
Pay 1 wed Ay d 
3 J 12 3 Jd 
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Order Input to MRA In MRA Transfer 
from MRA 
(continued) 
isu g3 - p. APs = p; = p 
3 


iy a ee “> ] + 


ime bl ae eae 


Q-0 


idv al -L 


its 25 





This group of orders will first convert the value of B currently 
stored in e6 from degrees into radians and store the result in e7. The 
last two orders will place the initial value of 862 which is assumed as 


unity, in g9. 
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Order Input to MRA Tn MRA Transfer 
from MRA 

ica e6 B B 

fie, el. Tt 1B 

; TB 

idv e2 180 80 

, ap mB 

ee 180 180 

1ca. Ing 1 at 

its ¢9 iL. 857 


This group of orders will calculate various parameters, such 
as cos B, sin B and store them temporarily for future use in the calcu- 
lation of end moments, The first instruction selects a third counter 
to be used for this and subsequent calculations, This is done for the 
Same reason as before, to permit this sequence to cycle three times, 
once for each inboard span, without affecting the cycle count of the 
previous loop. The second instruction sets the number of times the Loop 
is to cycle. This number equals the number of spans since m3 takes its 
valud from d3 which is the flad assigned for the number of spans, The 
manner in which the subroutines are utilized is shown here and their 


advantage becomes evident. 





m3, 


nd, 


Order 
ase he 
ale 
ica e/ 
isp s2 
deus ae 
ica e/ 
isp ls2 
tbs eZ 
ica ef 
isp s3 
its u3 
ica ef 
isp 5s3 
its ui 
lea 13 
isu wu2 
its u5 
Led Ul, 
rie WS 
2G “U6 


Input to MRA 


S:inh BS 


~sin By 


cos By 


% sinh By 
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In MRA Transfer 
from MRA 
By 
cos BS 
cos BS cos By 
By 
Sin By 
sin By sin PS 
BS 
sinh By 
sinh B. sinh BS 
By 
cosh BS 
cosh By cosh By 
Sinh BS 


sinh P,-sin Bo 


sinh 6,-sin By sinh B,-sin Bo 


cos BS 


cos By Sinh Bo 


cos By Sinh Py cos BS sinh Ps 





Os 


This group of orders will calculate the parameters sl, s2, s3, 


and store the results in g6, g7, and g& respectively. This operation will 


be repeated three times, once for eacn span, since the first two orders 


of the previous section are still in effect, 





Order Input to MRA In MRA Transfer 
from MRA 
ica uw2 sin By sin BS 
imr ul ~.cosh B, sin 8, cosh B, 
isu ud -cos B, sinh By sin B, cosh 8, -cos B, sinh B, 
Sin By cosh B,-cos 8, Sinh By 
idv ud + sinh B,-sin B SE B,-sin 5 
sin B, cosh B_-cos B_ sinh B 
Sueno gene. "4 
1 iD 
ica u3 Sinh BS sinh By 
isu wu2 -sin B, Sinh B,-sin B, 
imr e7 x By p, (sinh B,~sin B,) 
its ud B, (sinh B,~sin B, ) B, (sinh B,-sin B, ) 
ai ors gwen HE -CcOos By -cos Py 
imr w % cosh 8, -cos BS cosh BS 
jad lh9 +1 l=-cos Ba cosh By 


l=-cos Py cosh By 
idv ud +§,(sinh 8 


‘idan By) 6, (sinh 8, -sin B,) 





am 


Order Input to MRA In MRA Trensfer 
from MRA 


(continued) 


l=cos By cosh By 


alias B (sinh Bj-sinp,) = 


1 
ica u3 sinh Bo sinh BS 
isu ue -sin Bs sinh B,-sin By 
its u5 sinh B,-sin B, sinh B,-sin ae 
ica eld 2 2 
imr e7 * By 2B, 
imr w2 % Sin Ba 2B, sin BS 
imr u3 X sinh BS 2B. sin BS Sinh B, 


eB Sin B, sinh BS 
idv wd + Sink B-sin Bs "sinh f,-sin B, 
2B, Sin By Sinh By 


its 28 ner 
sinh By sin B, 3 


This group of orders will calculate the end moments and the 
angle of rotation at the inboard end of the shaft selected, The previous 


results stored in g5, namely iE and previously calculated values of 


"0 
8559 sl, s2, s3 are combined into the formulas shown in the flow chart 


to determine the values at the end of the second section, Mos This 
process is repeated three times to produce M3) 9 the moment we are 


looking for. 93), is also determined although it is not utilized. 





(2-6 


The method of saving registers that is used in this section should be 
noted. mm), is in g5. As each successive moment is calculated it too 
is put in g5. The preceeding moment is lost by this process, but it 
is of no importance since only the final end moment is required. 

The last six orders have a different function. The first five 
of these are utilized to calculate the value of B, for the next span 
which is used in the calculation of the moment and sl, s2, s3. The 
relation of B and B., is shown in the flow chart, 

The last order is used to cycle the program back to h8. ict 
causes the computer to return to h8 and at the same time increases the 
index of the cycle counter by one. When the index equals the number 
previously set in the criterion by an icr instruction ict will reset 


the index to +0 and pass control to the instruction immediately following, 


Order Input to MRA In MRA Transfer 
from MRA 
ica 95 Mo Moy 
its gl0 Moy m7 
= 7 _ 
ies g5 ms mM 
imr g6 X55 =m 85 
its ul ~My My 
ics f2+c = =C 
| 2} 
8 XS -C S 
a 3 a3 
imr g9 K 857 “C5 3° O1 





yaad 


its 


es 


idv 


isu 


its 


ica 


idv 


its 


ico 


g5 


gl0 


g7 


f2+e 


g6 


g? 


f/+e 
e6 


el 


e2 


e/ 


h8 


736 


Input to MRA In FRA Transfer 
from MRA 
(continued) 
ary -C Ts 
ee a 301 "O11 
- TONS 
Ye ol 1 1. 2 
aon econ! 
AY5 —My7°0 
aos 
<6 “O12 
52) C 
a 
Ee Mo1°0 3 
iL C al 
a 
Noisy a 2 
C 1 12 
a 
C C 
B B 
X C 
B BP 
x TY TIC 
BP 
‘ Tt 
£180 a C8 
Tl 
180 “pP B, 


This group of orders will print out the contents of g5, which 


are the moments at the inboard end of the shaft, 


The first instruction 





7h. 


is used to return the counting of cycles back to the 135. This is 
done since we want 135 end moments printed, or one for each 8 value 


used. 


Order Input to MRA In MRA Transfer 
from MRA 

isc l 

ica g5 M3), My), 


jmoatp 1234.56 % Ore 


This group of orders does several things which will be dis- 
cussed in the order of their occurance. 

The first three orders are used to increase the value of £8 
in e6 by 2° before the next cycle starts, This operation is performed 
135 times also, the B starts at 90° and goes to 360° in steps of 2 

The fourth order is used to break the normal sequence and 
add a correction before returning to the main program, a process known 
as patching, The correction added is given below, The first instruc- 
tion is merely a new flad used to get the number +1.0 in the computer, 
The second instruction increases the index of counter 1 by one and per- 
forms the same function as explained above, namely, permit this part of 
the program to cycle 135 times, The third and fourth instructions are 
merely to print a series of zeros in the print out to separate the 
results into blocks of 135 lines of freavency vs. moment for each order 
being investigated. The last instruction is to pass control to the 


next patch, 





Roe 


Order Input to MRA In MRA Transfer 
from MRA 


ica e6 B B 
iad e10 +2 B+2 
its e6 B42 B42 
isp 3h? 
isp hl2 
ict h6 
h9, i STOP 
START AT hl 
1h91+1,0 
ee in 
ica ih9 +1,0 41,0 
jimoat12 3liec 


isp-2h9 


This group of orders will do several things, First it passes 
control to the flad hl2 which introduces the subroutine of resetting 
the value of B to 90° at the completion of the program for one order. 
Then counter zero is selected to begin calculations for the next order, 
The instruction OUT is used to enter the Whirlwind mode and the next 
three instructions are used to increase the address of ||h6 by two. This 
is normally done automatically by the use of counters and #c at the end 
of each instruction. However, difficulty was encountered and this 
Manual" stepping was used to save time. The last two instructions are 
used to return to the CS mode and the main program at one instruction 


before h9, 





76. 


Order Input to MRA In MRA Transfer 
from MRA 
WL, +90,0 


~2n9/isp hl2 


hl2, ica hill +90.0 +90..0 
its e6 +90.0 +90 ,0 
asc) «0 
OUT 
calth6 285 285 
ade 11 +1. 286 
ade 11 +) 287 
tshh6 287 bh6 
IN 
isp h9-=1 


The main program is reentered as explained above and the last 
three instructions listed below are performed, The first increases the 
index of counter zero by one and when the entire program has been repeated 
three times passes control to the second instruction which stops the 
computer. 
The last instruction is used to start the calculations at the 
proper place. It is always the last instruction fed into the computer 
Since the computer scans the entire program, locates the flads and stores 
the information before it starts. Jf this instruction were anywhere 
else the computer would start computing before all the information was : 


fed into it and properly stored, 





ie 


Order Input to MRA In MRA Transfer 
from MRA 
ict h6 
ho. SLOP 


START AT hl 








78. 


PROGRAM OF STRESS CALCULATIONS OF NON-SYMMETRICALLY 
LOADED TRANSVERSE WEB FRAME 


The working program of the web frame stress calculations will 
consist of a series of sub-routines joined together by a sequence control 
program (see Figure V). To facilitate programming and program alteration, 
the instructions and input data were placed in four sevarate sections 
and these sections were represented by the flexowriter tape number that 
contained the instructions on data. Scantling data and other information 
relating to a particular problem were placed in tape 359-217-100, The 
sub-routines were placed in 359-217-101, and, the sequence control pro- 
grams were placed in 359-217-221 and 359-217-222. Calculations for a 
different though similar shin will require the commlete chanring of 


tape 359-217-100. Tape 359-217-101 is seneral for similar problems. Tapes 





359-217-221 and 222 determine the nature of the calculations, that is 
Whetier 2 check calculation or a désien calculation will be calcvlatec: 

In the instructions for the input data in tape 359-217-100 abso- 
lute addresses were assigned to flads. The absolute addresses were necessary 
Lom lock in,..bi out procedures used in the auxiliary storarve drum. The 
use of flads is still retained due to the ease of writine and caterorizing 
letters. The flad and absolute address assignments for storare of 
generated data were also placed in tape 359-217~101. 

“The following list of data are the specific input data for a 
DP type vessel. Due to the quantity «f numbers in each ategory, only a 


sample of the total group is presented below. 





32 | x1, +00.00! 
+3,266! 
+6.508! 
47 ALT 


(oe 


fe Tape 359~217-100 
(2,6) 


Assigns the first value of the x-coordinate of 
the shell to register 32 which is the first 
register in the case memory, Note that the 
coordinates are in feet measured from the origin 


at the keel. centerline, 


Total of 25 double length numbers which occupy 50 registers. 


+0 .000 
82 | yl, +00.90' 


25 numbers 


+2633 
332| v1, +2) .00 


25 numbers 

#12 .900 
282|b2, 409 3125 

25 numbers 


#0,2500 
482|d3, 412.7 


25 numbers 


ree. 
432\bh, +0.500 


25 numbers 


*0,250 
382|b5, +0.7500 


25 numbers 


40,375 


Value at the main deck centerline. 


First value of y-coordinates of the shell. 


ye-coordinate of main deck centerline, 
Value of web depth (a) at the origin measured 


in inches, First value stored in register 332. 
Value at main-deck centerline, 
Thickness of web (t,) measured in inches, 
Width of flange (w,) measured in inches, 

Ae 


Flange thickness (4) measured in inches, 


Shell thiclmess (t) in inches, 





peer e>, 970.0 
50 numbers 


~97h,.0 
932 e2, +25.517 


978 el, +8.0 


pal = 25 
pa2 = 18 
pa3 = 
pay = 
pad = 18 
paé = 18 


80. 


Weight and shear loads at each station. the 

weights are determined in averaging the weicht of 
vessel and shear for one foot of leneth. the length 
per foot times the number of feet of cirth between 


stations gives the weight load per station. 


Heisht of water pressure head. In this problem 
the desipn water head was the height of the deck 
edge. 


Longitudinal distance between web frames, The 
stresses computed were for the loads on one foot 

of logitudinal distance, Each web frame is assumed 
to support the loads imposed in the longitudinal 


dTSetance ~Deuneen™ Lr ames « 


Present parameter which represents the number of 


stations from the orisin deck centerline. 


Present parameter which represents the number of 


stations from the origin to the deck edge, 


Number of stations from the origin to the bottom 


of the starboard stanchion, 


Number of stations from the origin to the bottom 


of pori stanchion when measuring counter clockwise. 


Number of stations between bottom and top of the 


Storoocarc. stancni on. 


Numoer of stations between bottom and top of the 


port stanchion. 


The following instructions are the sequence control instructions 


and serve to utilize various sub-rceutines to perform the calculations of 


stresses for a web frame of one particular set of scantlings and load 








Gils, 


conditions, To alter the web frame and to produce a rational design, 
another set of instructions designated tape number 359-217-222 will be 


required. 


fo 359-217-221 (Name of tape) 
(2,6) 
(Selects number arrangement of double length registers. 


This arrangement utilizes 2) binary numbers for the 
mantissa and 6 binary digits for the exponent of the 


number). 

Pl Sp 262) gl forms sin 9, cos © and 4s 

1s “e2e=2 ge~e forms horizontal and vertical water loads 
Ho w 

isp gl0-c gl0-2 forms the total vertical loads by summing 
the values of the weights and vertical water 
loads. 

isp ses ¢3 forms sin ¢ and cos 4. 

isp gli Forms the reciprocal of the moment of inertia 
(1/Tn 4 ), of the section (A, ) and distance 
from the shell to the neutral axis (d,, A ie 

OUT The following orders are block cut instructions 
for transferring information from the case 
memory to the auxiliary drum. 

ca cf 20h8 (Address of first register of the drum 
storage block), 

si 967 Selects case output routine to the auxiliary 
drume 

ca cl2 800 (Number of registers necessary for storage). 

Bor e32 (Case memory address of the first register being 
blocked out). 

IN 

isp g6-2 Offsets the x and y coordinates of the neutral 


axis from the shell coordinates, 





isp el 
isp g3 
ica c6 
ts tll 
ies 03 
isp e7 
isp eg9 
isp g23-2 
isp g2h-1 
ica dl+50 
Scag rl 
imr 4@3 
aia ioe OS ' 
atts) a) cat 
OUT 
ca cld 
si 963 
ca cl7 
bi 332 
Carell 


see 


Forms sin @ and cos @, and As of neutral 
axis. 


Forms sin ¢/, and cos ¢ of neutral axis. 
+00 +0.0 
+0.0 
+0.0 
Forms m/I at each station due to applied loads, 


Forms integrals of M9 Pas Ans R and 5, into 


O 
a matrix (See Figure VI). 


Uses Cvout's method to solve the matrix, 


Alters V and H to include redundants V=dos 
H+D ys R+V and V4+sS, 


~m, malt 


0 


et m_/I 


m/T, 


Forms momente at each station, This moment 
includes the effect of the redundants 


2348 - Drum address of the first register of the 
drum storage. 


Selects "block-in" procedure for transferring data 
from the drum to the case memory. 


+150 ~ Number of registers being "blocked in," 
Blocks in d oe and Une 


Address of the first resister in the core memory 
that will receive the first number moved in from 
the drum. 


+2718 - Address of the first register of the drum 
storace. 





oo. 


si 963 Selects "block in" routine from the drum 
storage. 

ca cl18 #50 - Blocks in 50 registers of information 
of ay, 

iba 52 Address of the first register in the core memory 
that will contain day A 

IN 

p2, isp g25-1 Forms 01, J, T, and J” for all stations around 

the girth, 


The following instructions, called tape 359-217-222, perform the 
stress calculations for a web frame for a given set of scantlings and 
compare the computed stresses with allowable stresses, Jf the computed 
stresses do not fall within the range of allowable stresses, then the instruc- 
tions below will alter the scantlings of the web frame to produce the 
desired stresses. For this program the scantlings are altered by increasing 
or decreasing the width of the flange. Notice that tape 359-217-222 repeats 
all the instructions contained in 359-217-221 and in addition contains a 
method of altering the flange. While this may seem unnecessarily repetitious, 
two complete sequence programs were deemed more convenient to handle and less 


susceptible to error than one partial tape with a connective addition, 


fe 359-217-222 


(2,6) 
Ole Tsp ei gl forms sin 8, cos @ and As, 
isp g2-2 g2=2 forms ae and i 
isp glQ=2 glO-2 forms is 
isp 23 g3 forms sin ¢ and cos ¢, 
isp eh gl forms L/Ty ys A,, and Aya. 


OUT 





ca cB 
si 967 
ca cl2 
bo 32 
IN 
isp g6-2 
isp el 
isp g3 
ica c6 
its ti 
its b3 
isp e7 
isp @9 
isp g23-2 
isp g2l-1 
ica di+50 
its t) 
imr d3 
its b3 
isp g7 
OUT 
ca c10 
si 963 


Ga’ cl 


8h. 


2018 - Drum Address. 
selects drum output routine, 
800 - Nvmber of registers being transferred. 


Core memory, address of first register being 
blocked out, . 


Forms (xy 4s Yuya) 
sin 8, cos 0, and As of neutral axis. 
sin ¢ and cos ¢g or neutral axis, 
+0.0 +0.0 
+00 
+0,0 
Forms m_ /I at each station 


Forms integrals of m,, p., q., Rand S. (See 
Figure VI). 00 


Matrix solution. 


Altérs loads to include redundants. ds Po» R and 5, 


paul +m 0 
m, 
1/Ty M/Tp 
m/I, 


Forms moments ateach station, 


2348 - Drum address, 
Selects "block in" from drum internal routine. 


+150 ~ Number of registers being blocked in. 





| ops Ya 


“ers aaa oie 


si 963 
ca cl18 


isa {32 


isp g25-1 
OUT 


ca cl9 


si 963 
ca cl8 
ipsa Oe 
20 976 
ep p2-h 
spO 
IN 

Lea ue 3 


US US 


ieoeranle 


85. 


Address of first register in core memory to 
receive block in from drum, 


27),8-Drum address of first register in 
MHlock in" group. 


selects "block in" internal routine. 
+50-Number of registers being blocked in. 


Address of first register in core memory to 
receive blocked in data. 


Forms 2 a o. and co for all station. 


+21,:98-Drum address of first register to be 
blocked in. 


selects "block in" internal routine, 
+50-Number of registers being blocked in. 


Core address of first register receiving blocked 
in data. 


Forms cornter co limit number of alteration 
cycles to 3(1 + number in 976). 


Stops computer after completion of 3 cycles 


-1,0 


Note that the above counter was used for the test runs of this 


In any case the alteration procedure should be limited to pre- 


vent the possibility of the computer obtaining an unreasonable answer, As 


an example, should the web frame be to large the web flange can be reduced 


by predetermined amounts, However, large reductions could give a negative 


length of flange. 


A future programmer may want to change the value used 





86, 


in the counter and therefore the number of cycles. However a counter must 
be used in all cases. 

Now it becomes necessary to compare the values of the computed 
stresses With allowable values. To do this the computed values must be 
plus and a small program must be written to insure that all values are 
plus. Also, the maximum stress value at each point of equal scantlings 
must be determined. Remember that the web frame is symmetrical port and 
starboard, but the values of stress are not symmetrical. After determining 
the maximum value of stress at each station, the values are compared with 


allowable values and the scantlings altered, 


ict pal 
p2, ica e5+c (ce) 7 (starboard side) 
icp ml dle q is minus perform sub-routine ml, If 
plus, do next instruction. 
oS. ul +7_ 
cm 
ica e5+pb8+pb8 +7, +7. (port Side), 
ae mL. Checks sign of 7° 
6S Ue ry ga 
ae 
isu ti -7 (starboard) T.(port) - 7.(starboard) 
ae cmc Checks sign of 7,(port) - 7 ,(starboard) and if 
minus, goes to m2, 
ica t2 +7, 
“ £ +0, 
p3, isu ¢20 -(+20,000) T., ~(+20,000)(Lower value of stress 
La) 
cp m3 lng (7-20, 000) 16 Minus, po to m3 to reduce 


scantling values. 





isu c2l 


ED pl 


ica b3+c 
tad cl 


its b3+c 
aeeecs 
its t3 
pl, OUT 
ca p2e+3 
Suc] 


ts p2+3 


leu Hc 
2a, US 
icp p6 

OUT 
eaccly 
si 967 
ecucl'd 


bo 182 


ca cS 


Si 963 


€a~ ci? 


87. 


~8000 (J,-20..000)-8000 (8000 /in* 


is the range of stresses). 


If value is minus, perform sub-routine pl which 
moves forward to the next station, 


+W, + 
+25 W pt od (Increases flange width). 
+, 
We e 
-(-1.0) $1.0 
+10 


ica e5+pb8+pb8 = =ica e5+pb8+pb8 
-2 ica e5+pb8+pb8=2 


ica e5+pb8+pb8=2 


+2198—eDrum address register. 
Selects drum block out procedure. 
+50-Number of registers being blocked out. 


Core memory address of first register of group 
being blocked out, 


+20):8-Drum address register of first number in 
a group being blocked in, 


Selects drum block in procedure, 


800-Number of registers being blocked in. 





GG. 


oe) iol +800 +800 

ad cl7 +150 +950 

Ed sc2 +30 +980 

ts t8 +980 
ea. bo +980 

si 963 Selects drum block in procedure, 
CSCI +800 

bi 980 


The two block in procedures put numbers and instructions 
back into the core memory in preparation for another cycle of stress 
calculations of the altered scantlings, The block out instruction 
prior to the block in instructions puts new values of web flange 
width into the columns of numbers on the drum. The first block in 
instruction moves numbers into the core memory. The second block 
in instruction came as a result of "debugging" operations when it 
was determined necessary to replace all the addresses in the in- 
structions back to their starting value involving the port side it 
was necessary for the programmer to write instructions that would 
decrease the addresses of the registers containing the starboard side 
values. In other words it was necessary to move up a column of 


numbers in the core memory. At the beginning of the second cycle 





89. 


the instruction addresses were in error by a number equal to the number 
of cycles in that particular sub-routine. By blocking in the original 
instructions which remain on the drum into their assigned positions 
all the addresses were reset to their original value, 

For cycles within a basic cycle of operations it is necessary 
to write a reset routine that returns the addresses back to their 


original values at the end of each inner cycle. 


IN 

isp gh Starts cycle at she 

isp pi+l11 These two instructions skip the block out 
procedure in the main cycle. Since the block 
out has occured during the first cycle, the 
values remain in the drum and need not be 
blocked out again. 

ore wo TOr 


The following sub-routines are used in conjunction with 


the comparison and alteration sequence. 





90 


ml, ita ml+2 
imr c3 x-1 0 (-1.0x-@) = + ie 
(Makes J , plus). 

isp 

fie Jeans. +7, +(,(starboard) 

isp p3 

Moca, O5+C +Up +W 

ol CL = Wart 0d 

a 

its b3+c | Wp ed 

isp p3+7 

cl9, +2),98-Drum address. 

c20, +20,000.0-Lower limit of allowable stress 

value. 

ez +8,000.0-Range of allowable stress values, 
Note that for this problem a very wide range of values was 
used to prevent the possibility of the computer seeking 
an impossible solution, If the stress change caused by 
the flange alteration were to be larrer than the stress 
range, then a solution would be impossible with this program. 

976 =2 


i START AT pl 


Jn programming it is necessary to use various constants to change 


units, to form equations and to change generated values by arithmetic 





F1e 


manipulations. These constants must be assigned registers or flads to 
be of use in the program, The constants listed below were assigned flad 


cl, c2, etc., the c being useful in remembering the nature of the input. 


Tape fe 359-217-101 


(2,6) 

el, +0,0 

eee 430.0 A constant used to form the width of shell 
plating which is considered to be effective in 
forming the web beam, 

Coe 0 A constant useful in changing the sign of a 
value in the MRA, 

ch, +1.5 Used in calculation of the moment of inertia, 

c5, +61.0 Weicht of a cubic foot of sea water in 1bs/ft’. 

c6, +0,.0 Used to clear registers where the value in the 
register must be zero. 

Cia, +2 Whirlwind number used in changing addresses in the 
interpreted mode, 

c8, +2018 

ae bey A Whirlwind number used to assign addresses on 

3 «-4 7) 4 1 

ell, 427118 the auxiliary drum, 

cl2, +800 Number of registers in a blockout 

c13, +216 procedure. 

cll, +1 Whirlwind number used in changing addresses. 

15, *12,0 Converts feet or inches to inches or feet. 

cl6, +0.707 1/ 2 = sin 5° or cos lS”. 

cl7, +150 Number of registers in block out 

c18, +50 procedures. 


tl, through +8, +0.0 - Temporary registers for calculations. 





22 « 


The following preset parameters are part of the innu.s 7 71. 


do not vary with vessel type. Note that each preset parameter liste” below 


is generated from the present parameters in the variable input tare. 


pbl = pal-l 
pb2 = pal=-2 
pb3 — pa2~1 
pb! = pa3 +pa3 -2 


Number of segments in the haif cirth of a 
web frame, 


Tiro less than the number of stations on 4 half virth. 
Used in moment formation. 


Number of segments between the origin and deck edge. 


Number of registers from the orisin register to 
the bottom of the starboard stanchion. 


Also contained in Tape 359-217-101 are the storage assion= 


ments for generated numbersa 


He? | ail... 40.0 
DITTO 
18  +0.0 


162 | a2, +0,%0 
DITTO 
18  +40.0 


232 | a5, +0.0 
DITTC 
18  +0.0 


532 | eli, +0,0 
DITTO 
18  +0.0 


582 | e3, +020 
1D Bald 
18 +0.0 


832 | d3, +0.0 
DLETe 
18  +40.0 


Sine of the angle storace, Note that this 
instrvetion puts zeros in 25 double length 
registers. 


Cosine storage. 


Girth segment lengths storage 


Storage for vertical water loads. 


Storage for horizontal water loads, 


Storage for the inverse of the moment of inertia 


/ 
ui Ina. 





oo. 


882| d2, +0.0 
DITTO 
18  +40.0 
732] dl, +0.0 
DITTO 
18 +0,0 


782] sl, Library sub=-routine for finding square roots. 


The remainder of tape 359-217-101 consists of the various 
sub-routines that perform the functions listed in the flow diagram, 
Ficure V. A discussion of each sub-routine is placed at the beginnine of 
the list of instructions that comprise the sub-routine. 

The form of presentation of the following sub-routines was 
found to be an aid in the organization and understanding of the problem, 
Since continuity of thoucht is so important to the writing of a work- 
able program, the increased effort in writing the flor of information 
into the MRA, accumulation of data in the MRA and the values trans- 
ferred from the MRA was considered worthwhile, Other methods of 
organizing material could have been used; however, the authors found 
that this method was both effective and presented the ideas in a 


convenient form, 





SINE 


The following sub-routine forms the sin @ and cos 9 of the 


frame segments by the use of the basic equations listed below. 


Instruction 


gl, 


Ter 


VES 


isu 


its 


its 


ita gl+21 


pbl 


xl+24+c 


xl+c 


tL 


pik 


t2 


As. = Axe + Ay* 
n n n 
Ax = XX 3 
Ay, Se Nan 
sin @ = Ay/As 
cos 8 = Ax/As 


-Vn-s 





Transfer 
from MRA 


Input to MRA 


This instruction transfers the address in the 

control element to the address of the instruc- 
tion in gl+2]1 and adds one. (See glossary of 

instructions), 


Sets the cycle count to the number of segments 
on the side of the frame, This routine will be 


performed pbl times 


+X +X 
n n 
“x a 
n-1 eet n= 
Ax 
n 
Nee ele Ax® 
n n 
Ax* 
n 





INS TLUC LOM ai £1ad sol. 


oD. 


Instruction Input to MRA In MRA Transfer 
from MRA 
(continued) 

ica yltc+c 5 +7, 
isu. yl¢c ay i at 
ts. t4 by, 
imr +3 Ay Ay* 

n n 
iad t2 Nee er - bx¢ 

n n n 
oF sl sl is a sub-routine that extracts the square root 

of the number in the MRA and returns the square 
root of the number back into MRA, 

its a5d+¢c As a + ae 
ica tl +AX +Ax 

n n 
idv ad+c TAX, dx /ds_ 
its a2+#2+c (cos @) = Ax /As 

n 9 dene 6 
ica +3 +hy_, Ay, 
idv ad+c ras, Ay Js, 
its al+2+c (sin 8). S by /As., 
ict git2 Puts gi+#2 in the control element and starts the 
cycle over again. 

isp Address put into this instruction by the ita 
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The sub-routine that follows forms the water loads from the 
origin to the main deck, In this sub-routine the water loads are assumed 
to be uniformly distributed along the shell segment and the center of 
pressure is assumed to act at the mid-point of the segment, After formation 
of the total load on the segment, the load is divided by 2 and half of the 
load is placed at each end of the segment, The approximation of a con- 
tinuous beam under a trapezoidal shaped pressure load by a uniform load 
is only valid when the shell segments are small in size. As the loads 


are formed, the program converts the loads into horizontal and vertical 









components. 
2 A eee 
x 
FulrZ Cos = Sp 
wn sive LY 
"2 
S 
(Ful), ‘i 
z \ 
(= 64(H- ty) As, © 
as “Onn 
SAAT Sold 
2 
C)2-) 
AS bic UaOn Input to MRA In MRA Transfer 
from MRA 
ita g2+18 (See glossary of instructions). 
Ler - pos Sets cycle counter for number of segments in the 


starboard side from origin to deck edge, 


Castes: =Vitc a Be 


pes Mewes Tne] Inti *n * aS hoy 





Instruction 
im cl 

iad e2 

imr c5 

imr ad+tc 
a oC. 

its tl 

imr a24+2+c 
its elj#2+c 
iad el+c 
its elj+c 
ica tl 

imr al+2+c 
its e3+2+c 
jad e3+¢c 
its e3¢c 
Can Gee 


isp 
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Input to MRA In BRA Transfer 
from MRA 
(continued) 
0.5 ~(y_3)47,,)/2 = an 
+H(water head) (ft) iy 
3 
46h /ft 6),(H T aay) 
kAs (ft) éhAs(H=y__) =F 
0.5 B /2 (assumes simple support) 
F Bie 
x io 
(cos Q) v4] (cos O) ery /2 i Me) Fae 
Ce ea 
V V +V = V 
Wael wip 0 W 
(v), 
n 
OE a2) 
F 
W_. 
+(sin 8), (, sin @) = (Ho) 
(Hh) 
(HO) (H_) tC) = # 
n 
CH 


cycle back to g2. 


Address provided by ita instruction. 
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Formation of vertical forces is accomplished by the pnrorrams 
listed below. Vertical forces are the algebraic addition of buoyant forces 
which act in the upwerd direction and the weight forces which act down= 
ward. Upward forces are taken as positive. Note that a separate vrogram 
is required for each side. The cee forces are symmetrical and can be 
stored in 50 registers. The weight forces are non-symmetrical and the 
weight for each of the stations must be stored. In the program for 
forming the starboard vertical forces, the addresses of the uo moves down 
a column of registers of increasing addresses which is the manner of 
operation of the eljtc address. To move the other direction up a colum 
of numbers a small Whirlwind prograrn is required to redvce the addresses 


by 2 for each cycle, 


Instruction Input to MRA In MRA Transfer 
from MRA 
ita (See glossary of instructions). 
icr pa2 
gl0, ica elite +V +V 
W W 
jad ed+c Ww, #V_ Wl = Vi 
Lt 
its eS5+c Vv, 


End of starboard side program. 


ica elec +V +V 
W W 
iad e5+pb8+pb8 Wy Voli = es 
its e5+pb8+pb8 V 


t 


Goes out into Whirlwind to change addresses 
of weight loads. 


QUT 
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Instruction Input to BRA In MRA Transfer 
from MRA 
(continued) 
ca gl0+h iad e5+pb8+pb8 jad e5+#pb8+pb8 
su c7 -2 iad e54+pb8+pb8=2 
ts gl0+h jad e54#pb8+pb8~2 
td gl0+5 e54pb8+pb8-2 
IN 
ict gl0+2 
Sp 


After forming sin 9 and cos @ of each of the segments, it was 
necessary to find the sin and cos of the angle ¢ formed by the tangent to 
the curve at each station and the horizontal axis. cos g and sin ¢g were 
determined by averaringe the cosines and sines of the angles of the chords 
on either side of the station. While it was considered to be more accurat¢ 
to find the tan ¢ or the slope of the curve at each station, the possibility 
of tan ¢ = Ay/Ax becoming infinity rendered this impossible. Further, only 
sin ¢ and cos g were necessary for the purposes of the program and tan ¢ 
would have required additional storage and programming, This implies two 
conditions that do not exist. One, the ~ate of chance of angles should 
be equal and secondly, the ; 

(0, + @,) 


sin 3 


does not equal 
Sin 0, + sin @ 


ak 2 
2 
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and 
(0, + oy 


cos 
2 


does not equal 


cos 8, + cos 8, 


5 ° 
However, the resulting error will be small as long as the station points 
are judiciously selected to maintain small nearly equal angle changes, 
Should more accurate angles be necessary the programmer could utilize 
trigonometric identities and equal angle changes. 

An additional problem in this sub-routine was to insure that 
the centerline values of angles on deck and bottom represented 0° and 180° 
respectively. As was mentioned earlier the sin and cos of age were placed 


at the corner station. 


Instruction Input to MRA In MRA Transfer 
from MRA 


g3, ita g3+#20 (See plossary of terms). 


er pod 
Ga all4+24c #sin 9, +sin cal 
jad al+tl+c +sin 8, +sin 8, + sin 8, 
sin 8 + sin 8, 
amr ce! X0.5 a, sin ¢ 
its al+#2¢c sin ¢ 
ica a2#2+4+c #COS 8, +cos 85 


iad a2+h+c +#cos 8 cos 8), + cos 8, 
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Instruction Input to MRA In MRA Transfer 
from MRA 
(continued) 


cos 8 + cos 6 





imr cl. 0.5 5 = cos g 
its a2+2+c cos ¢ 
ict g3+2 


The following orders sets 0° and 180° in the | 
center line deck and keel and, also sets 135 
in the corner angle. 


ica c6 +0.0 +0.0 

its al +0.0 
its al+pbl+pbl +0.0 
26S C3 +1.0 +1.0 

its a2 +1,0 
ica c3 -1,0 ~1,0 

its a2+pbl+pbl -1.0 
ica cl6 +0./07 +0.707 

its al+pb3+pb3 +0.707 
ics cl6 ~0.707 -0.707 

its a2+pb3+pb3 =O. 707 
SO 


Since web frames are built up sections it is necessary to 
compute the moment of inertia at each station around the girth. The width 
of shell plating assumed to be effective in forming the web beam is 


thirty times the plate thickness. 





! ae 


em 1 


Cur 


dw 
1/2 ne |. MAL 


| 
ts 
See a be ee 


| 30ts —_—____—___~ 


The following sub-routine uses the relationships below in 





forming the area, moment of inertia and distance to the neutral axis 


of the web from the shell. 





ait a0 
A. “ Wot, mM. = (d + 5 )A I, = (g + ; )m.. 
Se a vee 
= Ss _— A = t—¥ 
oe Se TG SW “i 1.5 Ww 3 
2 e te 
— — —A eo =—— 
ae ey us 153" Os a ier 
- A +A 4A = = 
AL et ae : m, m.+m +m. I, Tat +1. 
m 
t, 2 
Va == Tha = FAL Sina 





Instruction 


eli, 
icr 


ica 


imr 


d 


d 


ica 


its 


ica 


iad 


its 


ica 


its 


ita g5+42h 


pal 
bl+c 


b24+c 


bl¢c 


cl 


blitc 
b3¢c 
t3 
bli+c 
cl 
bl+c 


t3 


b5+c 
b5+c 
c2 


t5 


Input to MRA 


oo 


xd, 


% 40.5 


+t 


XW 


+t 
X 0.5 


+d 
W 


+t 
xt 


% 30 


103, 


In MRA Transfer 
from MRA 
+d 
Ww 
aqatteA 
Ww WwW WwW 
A 
W 
4% 4, 
x = 
0.5 KA d_ m 
m 
W 
+t. 
ty x We = A. 
A 
+t, 
0.5 x te 
d_ + t,/2 
A(d_ + t,/2) = Hp 
Me 
+t 
s 
2 
(t.) 
2 
30(+,) =A. 
A 





Instruction 
imr b5+c 
lm - Gl 

its +6 
nea te 

iad ty 
1Su “tO 

ies. b/ 

ea bl 

iad t3 

iad +5 

its d2+c 
ea) Gs 

idv d2+4c 
its dl+c 
eo ica. t2 
imr bl+c 
idv ch 

ts. cl 


10h. 








Input _to MRA In MRA Transfer 
from MRA 
(continued) 
A x 
fee S 
ee 
X05 D = ae 
m 
S 
+m m 
W W 
+m, m + Mp 
“mn. m +m.-m « 
M 
15 
+A +A 
W WwW 
+A, A A 
+A A +A 4A <A 
S S WwW 16 
A 
+ 
+i sl 
M 
; 18 
A — oo qd 
os Nea 
qn A. 
+m +m 
W W 
A ae 
Xd md a 
W WoW 2 
we 
aS z= I 
I 





Instruction 


ica blitc 


iad? =bil+c 


jad tl 


ica t6 


idv ch 


imr b5+c 


lad tl 
ats | GL 
Tes. “Ge#c 
imr dl+c 
imr di+c 
iad tl 
ies tl. 


Tes C3 


Dey a aL. 


Its dc 


105. 











Input to MRA In MRA Transfer 
from MRA 
(éontinued) 
Ue Ue 
X05 t /2 
a. d +t,/2 
2 
X Mp m(d +.5t,) = A (da +t,/2) . 
+1 Tt, 
(I+t ) 
+m +m 
S S 
m4, ty 
F105 75 5 
wt 
Xt. 7 I. 
See tes 7 oe shell 
I nell 
“A, “A, 
Xdy A, “AL dy a. 
4 “AI a@ 
NA. NA. 
2 
OnE ete ee eee Gna. 
Iya. 
~(-1,0) +1,0 
sty oa. 1/Ty a 


Ty. 
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When working with beams, all forces and moments are considered 
to act along the neutral axis, To make the program more general, the 
original coordinate inputs were the x and y coordinates of the inner 
surface of hull shell. The coordinates of the neutral axis may be formed 
by offsets from the shell coordinates, The reader will notice that such 
a program allows alteration of the web section as desirec.to produce tne 


practical results, The offset distance is rar 


(Kiva, Yun) 
| 
MA. cosp!" On, 


Es 
dW.A. Sin g—< Xs, Ys) 


p 


it 


x - Oya. sin 4, 


Yuya. = Vg t hy, COS Be 


Ae 


The following program forms the coordinates of the neutrel axis. 
Since this program is a continuation of the above sub~routine, the counter 


instructions are still in effect. 


Instruction Input to MRA In MRA Transfer 
from MRA 

g6, ics dl+tc thy a. ~d, , (inch) 

idv cl5 +12 “dy , /12 (ft) 


43 ; es 
1ts tL dy ay [2 
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Instruction Input to MRA In MRA Transfer 
from MRA 
(continued) 
i al+c sin ¢ - ay 7 F 
Sa m shell 12 
: dy a,Sin # 
jad xlec +X shell we 1D = NR 
its xl+¢tc xy A. 
ics t1 o(-dy , /12 tay , 2 
: Neh, 
imr a2+c XCOS # nell 5" cos p 
any A 
—tiphe te 
Te Rosas Y shell Vs - 12 cos p = SNA. 
its yitc Yen 
ict gl 
Lep 


The moments at each station are formed by miltiplying the 
summation of the horizontal and vertical components of the loads by the 


vertical and horizontal projections of the sesment, 
ue SS 


—— eel 


(Xn ~Xo-) Vee 


i. 
(Yn 3 Yn— |} | 
Ca Yn-1) 


—— + — 





D>. Vn=r 
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» Vig (5%) i Hy g-Fp-y) ig ™, = ¢ 


The following sub-routine forms the moments at each station due 


to the weight and water pressure loads on the starboard side, Note that 


the first few instructions put +0.0 into tl, t2, t3. This was necessary 


due to their use in forming the summations, Jn addition when using this 


sub-routine to form ms tl, and b3 must be set to zero. When using this to 


compute the bending moment after solving for the redundant, the value of 


the redundant m, mst be put into tl, and m /19 in b3. 


Input to MRA 


Instruction 

gi, ita g8+35 

tem ppl 

ica c6 +0,.0 

its tl 

ats te 

its t3 

ica ed5+c ey 
n-2 

14a wt + V 

ats. tL 

ics xl+#2+c ~(x ) 

Moe 

n=-1 

imr tL x V 





In MRA Transfer 
from MRA 
+0.0 
+0.C 
+0.0 
+0,0 
ea 
n=. 
>." 
n= 


4, 
, 
= 


V 





Instruction 
Les: 03 
ica e3+c 
iad te 

Tes: ste 

ics yt2+c 
iad yltc 
imr t2 

Mea G3 

jad th 

ese Stl 

imr d3+2+c 
1ts b3+2+c 
ict g7+8 


109. 


Input to MRA In MRA Transfer 
from MRA 
(continued) 
(m_) 
+H +H 
W W 
n=2 n=1 
H H 
Ww W 
n=1 
>! 
Ww 
Tan H yan 
(7 y ane] ' Cares 2 Oy. nea] = 5 
n-1 n= 
x H A H 
} . W 5 W 
(m4 (m),, 
n=] 
(m_) (m_) 
(m_) 
x 1/1 m_/I 
mn /I 
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The following section of the sub-routine forms moments port 


side. Remember that the coordinates of neutral axis and the horizontal 


loads are stored for one side. 


The sub-routine mst move up the 


stored columm of numbers to form the moments, 


Instruction 


icr pal 


g8, ica e5tpal+pal+c +V 


jad tl 
its tl 
ica x1+pb2+4pb2 


isu xl+pbl+pbl 


its 


ica 


iad 


t3 


e3+pbl+pbl 


Input to MRA 


n=2 
V 


#x 
n+l 


-X 


n=l 


cal 


n=2 


n=2 


In MRA Transfer 
from MRA 
4V 
n 
n-1 n=] 
V4V oc V 
; n 
n-1 
V 
+X 1 
X41 *n eo 
n=1 
A V = Mm 
> n 
(m_) 
+H 
n=2 n=] 


HW = a 





Instruction Input to MRA In MRA Transfer 

from MRA 

(continued ) 

n-1 

res | te H 
n 
ica ylepb2+pb2 (475 Jay (49 A nel 
: - Y _ 
LSU yi+pbl+pbl (sn e é Ga EES (Yu a. | 
n=1 n=1 
ane a oA oe H = (may 
jad 13 (+m) (m + (mss =m 
n=1 n~l n 
iad ml+pbl+pbitc m aes Tl d= m 
n 
imr d3+pbl+pbl Le > m/I 
OUT The program now goes out into Whirlwind to 
change the addresses for each cycle. 
ca g8+3 ica xl+pb2+pb2 ica x1+pb2+pb2 
td 8+ Puts the address of g8+3 into the adcress of 
g84l; now reads isu xl+pb2+pb2. 

su c7 ~(+2W.W, ) ica x1--pp2+pb2=2 
ts g8+3 ica xl+pb2+pb2=2 
ca gS+7 ica e3+pbl+pbl ica e3+pbltpbl 
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Instruction Input to MRA In MRA Transfer 
from MRA 
(continued) 
su c7 ~(+2 WW.) ica e3+pbl+pbl-2 
ts g8+7 ica e3+pbl+pbl-2 
ca g8+10 ica yl+pb2+pb2 ica yl+pb2+pb2 
td g8+11 (Instruction in g8+1l is now isu yl+pb2+pb2 ). 
su c7 -(+2 W.W. ) ica yl+pb2+pb2=2 
ts g8+10 ica yl+pb2+pb2-2 
ca g8415 +imr d3+pbl+pbl +imr d3+pbl+pbl 
su c7 -—(+2 W.W,) +imr d3+pbl+pbl1-2 
ts 28415 imr d3+pbl+pb1~2 
IN 
£et (20 
isp 


Since all moments due to applied loads and forces had been 
formed, it was possible to perform the integrations of moments to ootain 
the slopes and deflections previously discussed in the procedure. At 


this point in the programming the redundants sPy9 and do must be con= 


0) 


sidered. will appear as a constant value of moment at each station 


M0 
throughout the girth, The moments due to Po and do depend upon the 


value of y and x at each station. Since the web frame was assumed to 


m7 .edundant _ Pon - Bowe 7 a 
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be symmetrical; the slope changes and deflections around the entire girth 
caused by the redundants will be double the slopes and deflections on one 
side, Fo: convenience in integration, the redundants will be integrated 
around the starboard side, 

The method of integration used was the trapezoidal rule which 
implies linear slope of the moment curve between stations on the moment 


curve, the area under the curve 


7 


(ay rn ee “Sn 
- ies I’n+1} 2 


n=1 


ds 
M My ds, Oh , | “2 
a= (2, + &, pe, + la we os 





As can be seen from the figure, the trapezoidal rule will cause small 
error providing the curve is smooth and stations are closely spaced, 
Despite the small error, the above method of integration was chosen 
because it is readily applicable to commer. programming and because 


random station spacing was necessary to approximate a curved figure 


oy a series of straight lines. 





1h. 


In the process of "debugging" the program several improvezents 
in the sequence of operations for the formaticn of the integrals were 
made. The results of these improvements and errors in addresses of the 
initial program for the integration procedures alters the procedure as 
written. The techniques used and explained in the following pages are 
correct. 

The major improvement was the inclusion of the integration 
routines under one sub-routine, All the integrals for My» Pos Ap and 
m_/I were performed under one cycle counter. Each of the stanchion 


redundant forces R and 5 were under a separate cycle. The routine 


yY - 
Po it: = 


was mot performed. To obtain the coefficient of Pos the coefficient of 


y 
fy ds 


was used. The group of instructions included under flac g? were added 


m, which is 


to insure that zeros were initially placed in the registers that were 
to hold the matrix. Also, the addresses for forming the port side 
integrals were found to be in error by two registers, 

The instructions listed below are correct and replace the 
integration routines listed under flads gl] throvch g2?2 of the detailed 
description. Unfortunately, these alterations were not made in sufficient 
time to allow the correction of the detailed explanation. These instruc- 


tions should be read from left to right, 





g9, ita g22 
ict g2+2 
imr a5+c 
imr xl+c 
imr xl+2+c 
jad dl+2) 
its tL 

iad tl 

eaip ice, tl 
imr yl+e+c 
its dl+l2 
imr cl 

imr b3+c 
iad t3 

its d1+52 
ica xl+2+c 
Lire ct. 

iad b3+pal+paltc 
its d1+50 
icayl+pb2+pb2 
imr cl 

imrb 3+pbl+pbl+c 
iad t1 

its d1+5) 
ts g18+2 

ts 219 

ts ¢19+10 
iG) ome 
icaxl+pbli+c 
mee U 

its d1+56 
imra5+pbli+c 
iad d1+56 


Les UL 


ACT. 3G 

Ler spol 

iad dl 

imr d3tc 

imr d3+2+c 
its di+2) 
ica yl+2+c 
imr a5+¢c 

imr yl+c 

iad tL 

gS aca bate 
jac d1+50 

its t3 

imr a5+c 

gl/, ica xl+c 
imr b3+2+c 
iad dl+5) 


imr a5+pb2+pb2 


ica c6 

gll, ica d3+c 
its dl 

its tl 

iad tL 

gl3, ica yitc 
imr d3+2+c 
iad dl+2 

its tL 

imr a5+c 

iad b34+2+c 
its d1+50 

ica ylt+2+c 

ah plage on B 

imr b3+c 

iad t3 

its di+5); 


vite sqm oa | 


£19, icayl+pbl+pbl imrb 3+pbl+pbd1+c 


imrb 3+pal+pal+c 


iad d1+52 

its tl 
imra5+pb2+pb2 
OUT 

ts g19+6 

su c/ 

Suc] 

ica dl+2 

isu xl+pbl; 
its tl 
icaxl+pbli+2+c 
imr cl 

its d1+56 
imryl+pbl+c 


iad tL 

its d1+52 
LOS +pb2+pbe 
amr 76) 

ca g18+2 

ts 219416 
ts g19+3 

ts g19+13 
its dl+10 
FuS. 6S 

imrb 3+pbli+c 
isuxl+pbh 
its t2 

ica tl 


iad di+16 


its ditc 

iad d3+2+c 
gl2, icaxitc 
ica xl+2+c 
imr a5+c 

imr d3+c 

its t2 

its d1l+2 

ica t2 

iad d1l+12 

imr a5¢c 

gl6, ica yl+tc 
imr b3+2+c 
jad dl+52 

oS 5 

imr a5+c 

¢18, icab3+pbl+pbl+c 
jad d1+50 

its tL 
imra5+pb2+pb2 
icsxl+pbl+pb1 
imrb 3+pal+pal tc 
iad d1+5u 

Su: cy 

ea £1943 

ca 219413 

IN 

g20,icr pad 
imra5+pbli+c 
iad d1+56 

its th 

imrb 3+pbli+2+c 
imrd3+pblitc 
its d1+16 


feay. te 
jad dl+16 
iad d1+26 
jad d1+26 
iad dl+36 
jad dl+36 
iad d1+6 
its d1+30 
its dl+3) 
its t3 
imrb 3+pb6+pb6 
isuxl+pb5 
its t2 
OUT 

ts g214+16 
its tl 
IGa uc 
iad di+18 
iad d1+28 
iad d1+28 
iad di+l8 
jad d1+)8 
iad d1+8 
its d1+ho 
its dl+)); 


imrd 3+pbli+2+c 
its d1+16 

its d1+26 

its d1+26 

its d1+36 

its d1+36 

its d1+6 

ica d1+16 
gels ier pao 
imra5+pb5+c 
jad d1+58 

its t 

imrb 3+pb6+pb6=2 
ca g21+16 

IN 
imryl1+pb5+c 
imrd3+pb5+2+c 
its d1+18 

its d1+28 

its d1+28 

its d1+)8 

its d14+h8 

its d1+8 

ica d1+18 
g22, isp0 


its 
Tea: 
1ca 
ica 


aes 


t2 
tl 
t2 
t3 
thy 


ev: Iie vl 
ict g204l 


its 


icaxl+pb5+c 


mr 


its 


imra5+pb5+c 


iad 


d1+32 


oak 
di+58 


d1+58 


ts g2l1+7 


ica 
iad 
its 
ics 
ics 
ica 
ica 
ica 
eu 


TtLy 


tL 
d1+18 
te 
tl 
te 
t3 
th 
al 
g2l+l 
d1+))2 


UG; 


imryl+pbli+2+c 
imrxl+pbl+c 
imrxl+pb)j+2+c 
imr tl 

imr t2 

jad t2 

ica d1+6 

ica d1+26 
isuxl+pb5 

its tL 
jicaxl+pb5+2+c 
imr cl 

its d1i+58 

su c/ 
imrd3+pb5+c 
its d1+18 
imryl+pb5+2+c 
imrxl+pb5+c 
imrxl+pb5+2+c 
abe gate 28 

imnr t2 

iad t2 

ica d1+8 

ica d1+28 





lt7. 


The results of each integration are stored to form the matrix 
as shown by Fig. VI, The following sub-routine performs the integration 


of the change of slope around the girth of the web frame due to redundant 





moment T° 
A A A 
0 "0 1 1 
of Ba nef pase am f G+ Bem fF +3 
0 O O 0 
Instruction Input to MRA In MRA Transfer 
from MRA 
gl1l, ita g15+9 (See glossary of instructions). 
Ter “pol (Sets cycle counter equal to the number of: 
stations on a side minus 1). 
ica d3¢c 41/14 +1/I 4 
iad d 
ia 342+c +1/1_ +1/T_ + 1/l 4 
imr ad+¢c x ds. (1/T_ + 1/1_,)4s., 


n=1 


ng) 
aad. di > Mt, + 1/1,_,)4s,_4 sy (1/I_, + 1/T,_, ds, 


n 
its dl (1/Z_ + 1/T,_,)ds_ 


The instructions below find the deflection in the y-direction 


due to the moment caused by the shear redundant dos 


6) 


: aa (x4) (x, ds 
Ay 2) x(q,x)ds = ‘. x(q,x)ds cs >). rs T oe 











Instruction 
gl2, ica xltc 
imr xl+c 
imr d3tc 
its tl 

ica xl+2+c 
imr xl+24c 
imr d3+2+c 
MG 1a 

imr a5+¢c 
jad d1+2); 
its dl+2) 


Input to MRA 


+X 
nwL 


n-1 


x Ii 


n=] 


118. 


In MRA Transfer 


from MRA 


2 
(x) /T 3 


2 
(x / 


x 
n 


x 
n 


eZ 
(x, )°/1, 
(x,)°/1, mi (x0) Tye 


(x )°/1, + Cae ds, 


n 
; com (x,/1,,)° uy (X49) Ty > Cx) pe) Tne os 


n 
a (5 )°/T (5) Tae Se 


The deflection in the x-direction due to m. was found by per- 


forming the integration 


0 


O ue on 
: aie of a Sadi 2D pea Eooa 9/8) im 
1 


The instructions starting with g13, perform this integration. 





Instruction 


gl3, 


imr 


its 


rca 


iad 


iad 


its 


of the moment due to Po = Poy around the girth. 


Q 
x= | 


ica yltc 


d3+c 


tl 


yl+2+c 
d34+2+c 
tl 


a5+c 


dl+2 


dl+2 


Input to MRA 


Tn-1 


x 1/l 4 


BAe 


x 1/1, 
S Vat n-1 


x ds. 


119 





In MRA Transfer 
from MRA 
*Tn-1 
Fyn T yA 
Vat one 
we 
ae 


(YS Tp /T7? 


Tne Tao 


Q= n 
S oy ene T3484 ; Oey LV n-1/ ea oe 
i I 


n 
F oF ee 
1 ds 


n 


The following section of the sub-routine performs the integration 


0 


0 





A 
f NG 4 ”n ES oor 
ds = 2p 4. ds = 2p + ) 
O Ai O/ I, 2 


eT 


alr 


/T a) 





Instruction 


gi, ica yl+tc 
imr d3+¢c 

its tl 

ica yl+2+c 


imr d3+2+c 


aes te 
jad tl 
imr abd+c 
iad dl+l0 
its dl+l0 


Input to MRA 


Tn-1 


x1/ oil 


bs ee 


x i/T,, 


*Vy/T 1 


x ds 


n-1 


120. 


In MRA Transfer 
from MRA 
Ty-1 
ey er 
SF sek 
ry 
y/t, 
y/1, 
Vel sae oa 
OE) net een 


n 
> Opal Teo! Tyee F859 ; Cz, LV y1/T 1) Sy 


n 
> el LSS eel )ds_ 


The x-deflection due to the axial redundant force Pp was 


Poy A 
(7e ds = 2Po 
O 


The instructions given in g15, performs the Po moment x-deflection. Note 


that the values of Yt and va! Tn formed in gllt, and stored in temporary 


registers is used in g15, thus saving operations, 





Instruction 


gi5, ica tl 
imr yltc 

its tL 

ica te 

imr yl+2+c 
ica tl 


imr ad+c 


jad dl+l2 


its di+12 


ict gll+2 


isp 


Input to MRA 


(Cycles back to g11+2). 


In MRA 


Vet nel 


ey) Mpa 


es iy I 
eae 
yt 


fa 2 
eae ey eee 


o2c% 


Transfer 
from MRA 


us 
Y ne? ee 


2 2 
2 Ee 


“> (y5/L_4¥°_4/Ty.4 48 
n 
; Z 2 
(YT ne! naw) n 
ab 


(Address set by ita instruction at 211), 


Due to symmetry of the web frame structure, the following re- 


dundant integrals are zero. Changes in slope due to Ios 


? 


o( 


OQ 
as { (4x) Cas +f 
A 





1228 
Vertical deflection due to Mm and Pos 


Q! ! A Qo! 


0 
ve Peay [i | Gees = { es 2) as a8 
A 0 


0! 


QB" (oy) A 
Ay = | mi ds = ro | “Bas + (-x)¥ds = Po QO. 
O A 


4 


Horizontal deflection due to dos 
oO! 


Ry (agx) . 
at T ds = "0 | Pe ‘i (-x)¥ds =) OF: 
0 A 


The sub-routines below perform the integration of moments caused 





by water pressure and weight loads, Note that two programs are necessary 
since 9 stations are used in the integration but only 25 x and y coordinates 
had been stored. The following program performs the integration of the 
moments along the starboard side to find the change in slope for one-haif 


the web, 


A 
m m 
n wee "n-1\ds 
se - ) Pes . F +z) 5 
6 n n=l 








Instruction 
gl6, ica b3+c 
jad b3+2+c 
im a5+c 

er. “cl 

jad di+50 

its d1+50 


of y m_/I to give the deflections in the 


half of the web. 


Instruction 
pije ea ye 
imr b2+c 

its t3 

ica yl+e+c 

imr b3+2+c 

Led 55 

imr a5+¢c 


Pa 


+m,_3/T 4 


n/n 


x ds, 


x 1/2 


n-1 


+ : (m_/T +m 
i 


In MRA 


+m 4/Ty 


(m_/T +m, 4/14) 


124 


transfer 
from MRA 


(m /T+m__,/1,_,)ds,, 


(m_/T +m, _,/1 448, /2 


n 
s/ T-)48,/ . : | (m,/ T/T )85,/ : 
ll 





iL 


n 
(m /I_+m__,/2,_,)ds_/2 


The segment of the sub-routine that follows forms the integration 


«ds 


Inout to MRA 


Wel 


xX mo 4/T 4 


we 


xX m/I 


Yy-1 m3/t 1 


In ‘MRA 


Wel 


Fina yea! pd 


2 


vn m/i, 


(ym /I 


ete? iy 


x-direction around the starboard 


Transfer 
from MRA 


r 
Yn-1 mo3/T 3 


n-1 m/T 1 


nn] m_;/1,_, ds, 





12 


Instruction Input to MRA Tr, MRA Transfer 
from MRA 
(continued) 
ime cl x 1/2 (y_, in pf tT ® ae ,/T__,)6s_/2 
n=] 
jad dl+52 ‘2 (y_m_/I at n-ne! nn 98 ,,/2 
ta 
1 
> (yn /L ty ann) 88y/2 
1 
n 
its dl+52 > (y_m_/L4¥ 3 1/11 )45, 
zh 


To obtain the deflections in the y-direction on the starboard 


side due to applied loads, the integral of x m_ /I was formed. 


A n 
poke a 1*n AL Mn oe 
= Pat! Al = eee ees i Sty 
ay = f qids - eee 
O if = me 
Instruction Input to MRA In MRA Transfer 
from_MRA 
Clo ica GAl+c +X +X 4 
imr b3te xm 3/T see one 
its +3 et 
ica xl+2+c +X x 
n n 
imr b3+2+c xm /i x m/l 
iad 3 +x, 7 y/t-3 (x *n-1 m/t nL * xm /t) 





Instruction 
imr ad+¢c 
ame ‘cl 

jad dal+5h 
its dl+5) 
ict gl6+2 


1256 


Input to HRA In MRA Transfer 
from MRA 
(continued) 
x ds. (x 4" 24/1 ae m /1 )ds | 
Roali2 a ee ae naa )as_/2 
n=1 
(x Te ytx,m Was )ds /2 
af 
n 
> (x ym /t tx 0 tn )ds ie 
1 


n 
> (et Ter Xt! tn) 98/2 
aE 


(Cycles back to gl6+2), 


To perform the integrations on the port side, the address changes 


if the x coordinates, y coordinates and girth segments ds must move up a 


column of registers, 


This is accomplished by subtracting 2 from the address 


for each cycle performed. 


Formation of the change in slope due to applied loads (port side) 








126s 


Instruction Input to MRA In MRA transfer 
from MRA 
icr pbl (Sets cycle counter to number of points in a side) 
g19, ica b3+paltpaltc +m__5/T 4 m-1/ 1-1 
jad b3+2tpalt+pal+c = +m_/T (m /T+m 5/17) 
imr a5+#pbl+pbl xds_ (m yt Ee _3/4,-1)4 
imr cl x1 /2 (m fi: +m _,/1, _, ds,/2 
254n-1 254n 
iad d3+50 > (m /I +m 5/Ta4 )as_/2 > (m_/ I +m /T,_,)ds_/2 
2o+n 
its d3+50 > (m tn +m Bea 04s, /2 


Formation of horizontal deflection port side 


25+n 


Ax = Pm a » Way n 
a i 2 
A 


Instruction Input to MRA In MRA Transfer 
from MRA 

g20, ica yl+pbl+pbl ea ere 

imr b3+paltpalte = xm, _/T 4 Ft nn! Tne 

its ee Ty-1 ra en 


ica yl+pbl+pbl-2 iss en +V 


Leis 


Instruction Input to MRA in MRA Transfer 
from MRA 
(continued) 
imr b3+pal+pal+2+c = x m,_/ 1, yn / L 
jad tl mV pe/t y ML nn ned 
imr a5+pbl+pbl xds_ (ym TAY n/t) 8, 
imr cl aye (ym /T AV t/t )98,/2 
254+n-1 
jad d1+52 (ym /T ty) ym _4/1,.1)48,/2 
as 
254n 
(ym /T AV pM / Ty )48,/2 


its dl+52 


Formation of vertical deflection port side due to applied 





loads. 
QO! Sen 
ca im He my? n on 
ty» | Base S (Sel, ony“ 
I T i 2 
A il ee n 
Instruction Input to MRA In MRA Transfer 
from MRA 
: + 
ica xl+pbl+pbl x1 +X 
imr b3+pbltpbl+c +m__,/1_4 Xe nnt 
ase AE? . 


Xa Td 





120. 


imseruct1on Input to MRA In MRA Transfer 
from MRA 
(continued) 
ica xl+pblt+pbl-2 +x +x 
+ 
imr b3+pal+tpal+2tc xm /1 xm /T 
i aE 
Ae eg n=l See a n’ n-1 ee ree 
imr a5+pbl+pbl+c xds_ (ym /T ty 1 /I_4) 
imr cl x72 (ym Val ae = ,)as,/2 
254n-1 
iad dl+5) (x m /T +x, ue ae = _) ds, /2 
a 
254n 


ae m _/ 1 +x yen ft, ds oe 


25+n 
its d145) Come Xe T ,)4s,/ 


The following program changes the addresses of the instructions 
in ¢l19 and g20 that utilize x, y and ds. The new addresses must move up 


a column of registers, 


OUT (Program goes into Whirlwind mode to change addresses), 





Instruction 


Ca 


Su 


ts 


ts 


ts 


Ca 


ts 


ts 


Ca 


ts 


su 


ts 


e19+2 
cy 
alk 
g20+6 
g20+17 
g20+3 
£20 

c/ 
g20+ 3 
22041) 
e2O0+11 
c7 


g20+1) 


IN 


Tew, 29 


isp 


Input to MRA 


imr a5+tpbl+pbl 


-2 WW. 


ica yl+pbl+pbl-~2 


-2 W.W, 


ica xl+pbl+pbl=2 


~2 W.W. 


129. 


In MRA Transfer 
from MRA 


imr a5+pbl+pb1 
imr a5+pbl+pb1-2 

imr a5+pbl+pb1=-2 

imr a5+#pbd]4nbl-2 

imr a5+pbl+pb1~2 
ica yl+pbl+pbl1~2 

ica yl+pbl+pbl~2 
ica yl+pbl+pb1-); 

ica yl+pbl+pbl-); 
ica xl+pbl+pbl-2 

ica xl+pbl+pb1-2 
ica xl+pbl+pbl-); 


ica xl+pbl+pbl-), 


(Puts computer back in interpreted mode). 


The above integrations provide the three equations for solving 


an unsupported ring frame of three degrees of indeterminancy. To solve the 


subject problem, the equations 





sa 





f —— = = —— 


3 





1306 


must be formed and solved simultaneously with the three previously formed. 
The designations Pie Bo By and 3h are the end points of the two stancnions, 
The sub-routine that follows determines the vertical deflection 


between B. and B. of the starboard 


1 le 
B B, 
A 7 hae 2S ete! A Mae nay? as 
ian — T i 2 
B B aig 
afk 1 
By | 
: miQo-x) as : m a1 (Xpa 7% Jas, 
a 2 I vs 
B n n-1 
a: 
Instruction Input to MRA In MRA Transfer 
from MRA 
ita g21+4)1 
ei aier pad 
i¢a xl+pbitec 2 | X-1 
isu. xl+pbl “x, (x, 77%,) 
its) oS (x) 
imr a5+pnli+c x ds (x, 47%, ds 
imr cl xe (x17, 48/2 
fies ced! Ci er leh | 2 
imr b3+pbli+c x (m/I) 5 (x, 47%, (m/T)_j ds, /2 ae 
n=2 n=1 


iad d1+56 7 > > a) (m/I)__,ds_/2 


i ai 





ies 


instruction Input _to MRA In MRA ‘Transfer 
from MRA 
(continued) 
n-1 
its d1+56 > (x __17%,,)(m/T)__,ds_/2 
aT 
ica xl+pbli+2 +x +x, 
isu xl+pbl = (x =x, ) 
acs? ty) (x =x ) 
eas 
imr a5+pbli+c x ds (x -x,)ds_ 
imr cl x 1/2 (x -x, )ds/2 
Les" (te ds/2(x_-x, ) 
imr b3+pbl+pb1+2+c x(m/ T) (m/ T) (x 7x, )ds/2 
B_+n-1 B. +n 


1 neha 

jad dl+56 +>" (m/I) (x,-x, )ds/2 >= Ga/D) (x -x,)a8/2 
B B 
is rh: 


B. #n 


al 
its d1+56 2 0/Malrayel 


The following instructions perform the integration 


as FE : _ 
y (x x.) ds = p ea x,) the Vey Xped x) os 
“0 i y i 2 i 3 
n n=) 


a 


Instruction 


ica 


imr 


its 


imr 


aac. 


its 


ica 


imr 


its 


imYr 


Taq 


its 


tL 
d3+pbl+c 
tee 


yl+pbl+c 


d1+16 


d1+16 


t2 


d3+pbli+2+c 


t2 


yl+pblj+2+c 


d1+16 


d1+16 


Input to MRA 


(x7 


x 1/ Beal 


x Vy 


B_+n=-2 
1 n 


n-L 


48/2 


132.6 


In MRA Transfer 
from MRA 
= 2 
Cans x, )ds_/ 
(x__47*,/1,,-1)98,/2 


(x,_17%,/T,_7 )4s,,/2 


Fe 1%) /7 a 88,/? 


si > Vn Mya %q)/ T198,/ > Vy © *-1" x, /I,, : 195,,/ : 
B 


a 


(x =x, )ds_/2 


By+n-1 
B 
1 
B,+n-1 
Yn &% = fee ak ds_/2 
aT 


(x -x_)ds_/2 


x(1/T), (x, -x_)/T_ ds, /2 
(x -x_) /1,ds /2 
xy. (x =x dy /Tas_/2 
Bj +n-1 Bi+n 
+. > _1(, —x a) Ty 1 ds Af? ¥en, tds /2 
B B 
1 1 
Bi+n 
> ¥ (x -x,) /T ds /2 


1 


1S Bi 


The vertical deflection between the ends of the starboard stanchion 


due to the redundant Vo force was found by 





*n- a n=-1 -x,) oe : Yn %y>%,) or 
I a 1 > 
n=1 
itis bruce bon Input to MRA In MRA Transfer 
from MRA 
meas = td. (x _47%,/4,-1/d8,,/2 (x__37*,)/0,_198,,/2 
imr xl+pbli+c x X 4 x21 '% =o ,/1, see 
Bi+n-2 B_+n=1 
i t 
jad dl+26 x a1 -x ft, _138,,/2 a _ {x -x tn _148_/2 
BS BL 
B +n-1 
1 
its d1+26 | x7 (X -x /t, _738,,/2 
a 
i - ees - uf 
ica t2 (x, x, )/ s /2 (x, x,)/ nas, /2 
: = - 
imr xl+pbl+2+c x xX, n'y x,)/1,as_/2 
pee 
yen al: Bi +n 
jad dl+26 : oe Pic. -x _)/4., _, 98/2 > x(x xy /E ds _/2 
2 oa 
Ben 
its d1+26 > x (x -x )/I ds /2 
nna noon 
B 


1 





134. 


The vertical deflection between the ends of the starboard stanchion 


due to the redundant colwm axial force R was found by 


=x) @ (x + d 
ir tee « Par sa Soa” Sa 








Instruction Input to MRA In MRA Transfer 
from MRA 

eae Go (x, 47%) (x, 47%,) 

: 2 

imr +1 X (x, 47x, )/T ds, /2 (x ox, )°/1_, 48/2 

Bj +n=2 Bj +n=1 
. ee 2 2 
iad d1+36 + (x 4 x.) / cee 2 > (x 47%, /T2198,/2 
a a 

Bo+n 

its dl+36 

its 3 > (x a ta) Pn ,as_/2 
ei 

ica th (x -x,) (x ~x,) 

; es aaa 

imr +2 x (x, -x, )/Tds_/2 (x, x.) /Tas_/2 

Botn-l B. +n=1 
2 
iad d1+36 re > (x, 47%.) /T__, 38/5 : (x -x Jy /T ds A? 
B 
al 

B yen 


its dl+36 xe “x, Ve /I, ds oe 





dS y 


The vertical deflection between the two stanchion ends starboard 


Side due to m,) Was performed below. 





2 (x-x, ) me cu (x -x,) ds . (x Eee ds 
04 fs T ieee 0), eee ame 2 i 2 
B n n-1 
1 
Instruction Input to MRA In MRA Transfer 
from MRA 
ica +1 (x__47%,)/1__4 48/2 (x__,-x,)/1,_498,/2 
jad t2 + (xox )/1,48,/2 (x__5-*,)/1,-798,/2 
+(x =x, )/Tds_/2 
B +n-1 | B +n 
aa hele +) | re ated oe a a 
25) ( au 
+(x - d - d 
‘i x /1, s /2 +(x x1, s_/2 
Bi+n 
its dl1+6 |pea¥a)/ Foe 
B 
+(x =x.) / | ds_/2 
26% e214] (Cycles back through the stanchion integration 


program). 


isp 


The following sub-routine performs the integrations around the 
port side which will give the total vertical deflection between the ends 
of the port stanchion, The entire sum should equal zero, since the stanchion 


is assumed incompressible. 


uses 


The instructions that follow perform the integration 


B B 


4m (xex,) 
oe a as - Es (am,) +p To 2 3) |S 


Ne” 


°3 3 
Instruction Input to MRA meals Transfer 
from FRA 
ita g22+18 
g22, icr pa6 
ica xl+pb5+c | a 
isu. xl+pb5 ee ee 
ts? %3 . (x 47%) 
imr a5+pb5+c x ds, (x, 17%, 4s, 
imr cl M72 (x 17% )ds_/2 
tone bil (x 47%, ds,,/? 
imr b3+pb6+pb6 x (m/I) 4 (m/T) (x47, /ds8_,/2 
B +n=2 Byen-L 
iad di+58 Sys Ggas%)86f2 > a ae ee, ds nf? 
°3 *53 
Bl+n=1 
its d1+58 | new, 723 ds e2 
a 
ica xl+pb5+2+4c , +X, +x 


isu xl+pb5 “x, x 7X, 





137. 


HnStruc Lon Input to MRA In MRA Transfer 
from MRA 
(continued) 
its th 2 
imr a5+pb5+c Kds_ (x -x, ds, 
imr cl Kole (x -x,)ds_/2 
its t2 (x =x, )ds_/2 
imr b3+pb6+pb6-2 X (m/T) (m/I) (x =x ds, /2 
B_+n=2 B +n 
iad d1+58 (m/T) (7 (x,7%, 48/2 (m/I) (x =x, )ds,/2 
B B 
3) 5) 

Bin 
its d1+58 (m/T) (x =x, )ds_/2 

B 

3 
OUT 
ca g22+16 imr b3+pb6+pb6-2 imr b3+pb6+pb6—2 
ts g22+7 imr b3+pb6+pb6~2 
su c7 imr b3+pb6+pb6-) 
ts g22+16 imr b3+pb64+pb6-), 
IN 


The following section of the sub-routine performs the integral 


B 
Ly (xe, ) Cee eu Geeseoaia 
ie). a = by prea rena te Holgi |S 


n=1 n 


Bs 


Instruction 


ea. * ti. 
imr d3+pb5+c 
oS, Ge 


imr yl+po5+c 


jad d1+18 
its d1+18 
eas ae 


imr d3+pb5+2+c 


its te 


imr yl+pb5+2+c 


iad dil+18 


its 4di+18 


y_, x 7%, )/ Ids / 2 


Input to MRA In MRA -rar.ster 
fror MRA 
(x, 17%, /08,,/2 (x 17%, )ds [2 
A/T) ry erate 
ey aes = i" 5/2 
XVn-1 ee Syne enn 
B_+n-2 B+n-1 


B B 
3 3 
Bo+n-1 
a ST ane 
: ¥ nat Xpar %)/T,_4d8_/2 
oe 
a 
(x -x%, Jds,/2 (x, 7%, ds_/2 
jaa te SA - ‘T as 
1/ n (x %,)/ eo wf 
ne aaa eee fe 
ee “D ae 
Vy Ty b%,7%,)/1,,98,,/2 
BL+n=1 Bo+n 


~ acre _-7 yr 3 /< 
Fray Xo 4/2198, /2 ean ie ae n : 


B. B 


Ng 





139.6 


The following section of the sub-routine performs the intecration 


x(x-x,) x(x ee) is Selene — ds 
~[ “PS Seed Poet as 


Instruction Input to MRA In MRA Transfe> 
a from gra 
ieee Seti (x, 47%, )/1, 7 98,,/2 oes “x /L 4 ds Al? 
j - I Se 
imr xl+pb5+c X Xx 4 a | x )/I__4¢ oe 
BL+n-2 B en-1 
iad d1+28 x (x =x )/I_ds/2 x nel -x )/1 a1 05 s_/ 
B B 
3 3 

Blen-1 
; ; = )/E, (as /9 
its d1+28 x {x pep te Sf 

B 
3 
~x,)/I d =x, )/l_d 
ica t2 (x, x )/ - s_/2 (x x ys 8/2 
imr xl+pb5+2+4c X xX, x (x -%, )/Tds /2 
Bl+n-1 B+n-1 
iad d1+28 x1 7%, / 1 ds /2 > x-1 Fp, =z, /1 ncs,,/2 
B iB; 
3 3 
eee 


p ; ; t ae 
its dl+28 x(x, x )/ nis, /2 


Bo 
2», 





TO 


The following section of the sub-rovutine forms the vertical 


deflection caused by the axial stanchion force redundant 5, 


B 


Gee) . 
Ay = sf pai cae ae eS 
amt 
3 eS 
Instruction Input to MRA 
iea-' v5 X17 
am y tal A ea ae 
B_+n=2 
Dee 
iad dl+))8 (x, 6 x) /I,_,98 /2 
ae 
its d1+)8 
ica th (x -% ) 
imr +2 X (x -x, )/Tds_/2 
Bi¢n=1 
jad di+l8 (x 7%) /T__ 48/2 
ae 
its di+)8 


B 


I 
n 


2 2 
= ny ee 


z 


n=l 


In MRA Transfer 
from MRA 
(x4 7%,) 


(x, 47%) Tg 38,/2 


+n-] 


(oe 


( x 7%) 


1%) Ty 88,/2 


Bl+n=1 


(x 4°%,)/ Tn1°8y/ 


= 


(x -x, )°/1,ds,,/2 


Bl+n 


: (x, _47%,)°/t,a8_/2 


B 


5. 


B_+n 


Dea, 


B 
: 





Bors 
tip » 


The following section of the sub-routine finds the vertical 


deflection caused by the redundant moment m,.. 


Instruction 
eae ty. 

wad. pte 

jad di+8 
its di+8 
ict g224L 


O 


B, By 
(x35 ) a (x -%,) Oa Se ) — 
ii cgi it ee Tay 9 
3 3 








Input to MRA In MRA Transfer 
from MRA 

(x__5-%,)/1,-198,/2 (x 47%, )/T,-198,/2 

(x, -x, )/1,ds_/2 (447%, )/T,27 48, /2 


+ (x,-x, /1 ds, /2 


B +n=-1 Bl+¢n 


a | 1%) STeaya)/y 28/2 
B 


B 
: +(x,-%,)/1, |48,/2 ° +(x x, )/Las_/2| 


B +n 
) 


Sz | a%, Ee 


5) » 9 
+x x/1, | ds__/ 2 


(Cycles back to 22241). 


The redundant stanchion forces affect the deflections for the 


web as a whole in much the same manner as the three basic redundant forces 


affect the deflecticn around the two stanchion sections, An analysis of 


Figure VI shows thet the integrals of Por 4 


0? and 5 for the traverse 


around the stanchion loop is the same as the integrals of R and S about 





LUC 


the entire web girth. Therefore, the coefficients of the redundants which 
are the integrals may be used in several places. As an exam-le the coeffi~ 
cient of the deflection magnitude due to My between 7 and Bs is the same 
as the coefficient of the angular rotation magnitude due to the redundant 


R about the entire girth. Moment due to R is R(x-x, ) and the rotation 


due to the moment is 











B B Qo! 
t 
: . (x-x,, a . (x-x, ) © (xex,) (x-x, ) 
i 7 eae = as + i as + i QS. 
O O B By, 


The deflection between BL and 25 due to ™m, is 


(x-x, ) 


9 |r 8s 


The following instructions utilize the computed coefficients to 


fill similar valuess 


Instruction Input to MRA In MRA Transfer 
from MRA 
B i 
ve 3 : 
ica dl¢6 (x-x, )/Tds +f (x-x, )/Ids 
B B 
a a 
B 
3 id ¥ 
its di+30 f (xex, )/Ide 
I 
B 
2 Po 
ica di+l6 . y(x-x, )/Tds - J y (x-x, )/Ids 


A) 


ee) 





Instruction 


its 


ica 


its 


nea 


its 


Loa 


its 


ica 


d1+ 32 


d1+#26 


d+ 3) 


a1+8 


a1+h0 


d14+18 


dl+)2 


d1+28 


Input to MRA In MRA 


(continued) 


B 


2 
wf x(x-x, )/Ids 


aL 


B 


h 
7: (x-x, )/Tds 


3 


1, 
f y (x-x, )/Ids 


Bs 


B 


mM 
f x(x-x, )/Tds 


* 


13. 


Transfer 
from MRA 


B 


2 
ay y(x-x, )/Ids 


a 


Lo 
f x(x-x, )/Ids 


B 
1 


yy 
(x-x, )/Tds 


es 


B 


) 
Be y (x-x, )/Ids 


ae 





Instruction 


its 


ica 


its 


eS 


10S 


ica 


Les 


es 


d1+ 32 


d1+26 


d+ 3) 


d1+8 


di¢ho 


d1+18 


d1+)2 


G1+28 


Input to MRA In MRA 
(continued) 
B 
2 
fo x(x-x, )/Ids 
By 
eh 
if (x-x, )/Ids 
3 
v 
f y(x-x, )/Ids 
= 
sy 
7 x(x-x, )/Ids 
B 


3 


13. 


Transfer 
from MRA 


B 


2 
i y(x-x, )/Tds 


al 


B 


2 
3 x(x-x, )/Ids 


B 
it 


hh 
(x-x, )/Ids 


ae 


B 


l; 
f y (x-x, )/Tds 


eo 





Wis. 


Instruction Input to MRA In MRA Trensfer 
from MRA 
(continued) 
B 
h 

its dle) x(x-x, )/Tds 

5 
isp 


The following sub-rovtine reduces the mtrix of deflection 
equations by Crout's method and finds values for all redundants, The sub- 
routine, designated by 52, is a routine contained in the library of the 
computer laboratory. The library routines have been tested and tried 
and are known to be correct. For this reason only the method of entering 


the routine will be given. 


Instruction Input to MRA In FRA Transfer 
from MRA 
ita g2345 
OUT 


g23, SP 82 
dl Address of first colum, 
Order of matrix. 
aN Number of columms of constants 
IW 

isp 

The instructions below alter the vertical end horizontal forces 
to include the recundant forces Pos 4 


O 


is applied as a load at its point of action. 


, @ and 5, Note that esch redundant 


pote = V eV - 
ey meee aes ope 





instruction 


ita 


gh, 


isu 


its 


ica 


isu 


its 


LJ 
C) 
9 


isu 


its 


1¢2 


iad 


g2h425 


ica e5 
di+5) 


— 
€) 


So 


d1+52 


e5+pbly 


d1+56 


e5+pbl 

e5+pb5+pa5+pa5 
d1+56 

e54+pbli+pa5+pa5 

e54pb6+pb6 
d1+58 

Bee neanG 


ab6+pb6=pa6—-pad 
d14+58 


pb&+pb6~pa6-pa6 


Input to MRA In MRA 
+V +V 
n n 


Yo “a4 


The following instructions form H-Do 


+H +i] 
W WW 


“Po H Po 


The following instructions form V-R 


+V +V 
n n 
-R V -R 
n 
+V 4V 
n n 
ree V +R 
1g) 
+V +V 
n n 
+5 Vo#5 
n 
+V +V 
ras 4V45S 


Hee 


eco 


Lrenste= 
from MRA 


Ve- 
n 4p 


H =n 


w  O 


al 


Ves 





thé. 


The Whirlwind instructions that follow change the adcress of 


geh+16 and g2)+18,. 


Instruction 


OUT 
ca g2h+16 
Su (Ces ~2 
ts g2l+16 
td g2h418 

IN 


ict g2h+16 


Input to MRA 


ica e5+pb6+pb6-2 


In MRA 


Transfer 
from MRA 


ica e5+pb6+pb6-2 


ica e5+pb6+pb6-) 


ica e54pb6+pb6-—); 


e54+pbo64+pb6—!: 


The following sub-routine finds the stresses at each station due 


to the axial loads (p), the shear loads (q) and the bending moments (m). 


After the various stresses have been formed the direct stresses and the 


maximum compression and tension bending stresses are combined to form the 


total stresses at the flance and at the shell surface, 


listed below are utilized in the sub-program,. 


p= 

q= 
Direct stress J. = 
Shear stress Tease 


Bending Stresses = 


YQ el yal 


H cos 6 = V sin ¢ 


H sin gd - V cos ¢ 


p/A 
q/A 
my /L 
shell, = m/I (dy, , +t.) 
ar = m/ - 
flange = m/I (tp+d_ 
oS 
Fa Sn 
7 OEE 


lhe relationsa.s 
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The following instructions use the adove relationships for the 


starboard side of the frame. 


Instruction Input to MRA In MRA ‘ Transfer 
from MRA 
ita 227 
oOo ica 1cO +0.0 +0.0 
Pos. Go 40,0 in “Gb 
Riess “op . +0.0 
acr ‘pal 
n-1 n=] 
a a 
es se . V _ V 
wi 
iad e5+c +V toe V 
n 
abs oC > V 
14) 


imr altc X sin g > Vsin 


n 
its tl > V sin g 
n=-1 n=} 
Tea (iu) + > H + H 
n 
jad e3+c +H + > H 
n 


US Gy: H 


i 


imr a2+c x cos ¢ H cos ¢ 


isu tl -V sin g > Hoos f- > Vsingep 





Instruction Inout to MRA In MRA Srace Tor 
fror MRA 
{eontinued) 

idv d2¢c arm p/A 

imr el X 8 8 p/A 

its xl+tc p/A = o- 
n a 

ica t8 + V fA V 


imr a2ec x cos ¢ > V cos g 


n 
its “td. : V cos ¢ 
n n 
ies £7 + > H + > H 
n 
imr al X sin ¢ H sin ¢g 
n n n 
tou td. - > V cos ¢ > Hsin 6 -pV cos faa 
idv d2+c ae q/A 
imr el X 8 8 g/A 
ats sic 8 afA = 7.” 


Note that, was formed using the entire areca as effective in 
Shears. This is an error which should be correcte? sirce only the area 
of the web is effective in shear. 


ics dl+c “dy A. LA. 





Anstruction 


isu 


iad 


its 


1ca 


aire al 


Jad 


its 


ae 


port side, 


Input _to MRA 
(continued) 

b5+c ae 
e15 X 12 
b3+c xX nf/I 
il X 8 
xl+c + 

a 
aS+c 
bl+c +d 
blitc +t, 
ee wdy oa, 
e1s es 
b3+c X m/I 
el X 8 
xil+c + 

a 
e5+c 
g25+3h 





Tn MRA = caneter 
from MRA 

—(dyy 4g +t,) 

“12(dy 4, +t.) 


“12 m/Idy 4 +t.) =, 


-96 m/I(d Gm 


+t ) = 
Ss sm 


NA. 


Oy Q) 20, 


Cd ttprdy 9 @ ) 


12(d ttendy y ) 


12 m/Td + temday. 9) 


96 n/Ud+tend , ) ~ 


CL vg 


fm “a f 


T + 


The following instructions calculate the stresses around the 





instruction 


Ov 


ge6, 


140 


aus 


inr 


ts 


eS 


isu 


us 


iad 


LAV 


its 


pb1l 


ica t8 


e5+paltpaltc 


t8 


al+pbo2+pb2 


ia 


t7 


e34+pb2+pb2 


t7 


a2+pb2+pb2 


soil 


d2+pb2+pb2 


xl+paitpaltc 


input to MRA 


n= 


+> Vv 


+V 


X sin ¢ 


n=1 


+H 


X cos p 


ee a 


gO 


In MRA Transfer 
from MRA 
n= 
+> Vv 





Instruction 


ea 


amy 


its 


e434 


isu 


ake aig 


ts 


ics 


isu 


ime 


isu 


t8 


a2+pb2+pb2 


tL 


tf 


al+pb2+pb2 


tL 


d2+pb2+pb2 
sl+pal+paltc 
dl+pb2+pb2 
b5+pb2+pb2 


b3+pal+pal+c 


xl+pal+pal+c 


a5+pal+pal+c 


Input to MRA 


(continued) 


nN 


> V 


X cos ¢ 


“Ci A. 
a 


Xm/I 


In MRA Transfer 
n 
DY 
ne 
cos ¢ V 


“dyn. 


= j 
(dy tt.) 


oa m/I(d., , +t) aes 0%. 


2 ae 
a sil S 


Ts 





152, 


ines traction Input to MRA In MRA Transfer 
rom MRA 
(continued) 
ics dl+pb2+pb2 “dey “dy y 
jad blj+pb2+pb2 +t, “dyn the 
jad bl+pb2+pb2 ede (dtterd  ) 
imr b3+pal+palt+c X mn/I m/T(drte-dy , ) = Orn 
iad xl+pal+pal+c 205 Gea, -On 
its e5+pal+pal+c ~~ 
OUT (The following instructions chanzve the addresses of 


the above program). 


ca g26+3 imr al+pb2+pb2 imr al+pb2+pb2 

su c7 ~2 imr al+pb2+pb2~2 

ts 226+3 imr al+pb2+pb2--2 
ts g26+17 

ca g26+7 iad e3+pb2+pb2 iad e3+pb2+pb2 

Sue 7c i, —2 iad e3+pb2+pb2-2 

ts g26+7 jad e3+pb2+pb2=2 
ca g26+4+9 imr a2+pb2+pb2 imr a2+pb2+pb2 

Siu Cy ~2 imr a2+po2e+pb2=-2 

ts g26+9 imr a2+po24+pb2—2 
ts g264l1h imr a2+pb24+nb2=2 
ca g26+11 idv d2+pb2+pb2 idv d2+pb2+pb2 

Su. C7 -2 idv d2+pb2+pb2-2 


ts 226411 idv d2+pb2+pb2-2 





oes 


Instruction Input to HRA Tn MRA Transfer 


(continued) 

ts g26+18 idv d2+pb2+pnb2—2 
ca g26+20 ics dl+pb2+pb2 ics dl+pb2+pb2 
su c/ -2 ics dl+pb2+pb2-2 
ts g26+20 ics dl+pb24+pb2~2 
ts g26+25 ics dl+pb2+pb2=2 
ca g26+21 isu b5+pb2+pb2 isu b5+po2+4pb2 
Su C7 -2 isu b5+pb2+pb2=2 
ts g26+21 isu b5+pb2+pb2=2 
ca g26+26 iad bl+pb2+pb2 iad bli+pb2+pb2 
su c7 ~-2 iad bll+pb2+pb2—2 
ts g26+26 iad blj+pb2+po02~2 
ca g26+27 jad bl+pb2+pb2 iad bl+pb2+pb2 
Su, "C7 -2 iad bi+pb2+pb2=2 
ts 226427 jad bl+epb2+pb2=2 

IN 
ey ice6 


gc?, isp 





Aopendix C 
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